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Abstract 



We study the moments of equilibrium measures for iterated function systems 
(IFSs) and draw connections to operator theory. Our main object of study is the 
infinite matrix which encodes all the moment data of a Borel measure on M.'^ or C. 
To encode the salient features of a given IFS into precise moment data, we establish 
an interdependence between IFS equilibrium measures, the encoding of the sequence 
of moments of these measures into operators, and a new correspondence between 
the IFS moments and this family of operators in Hilbert space. For a given IFS, 
our aim is to establish a functorial correspondence in such a way that the geometric 
transformations of the IFS turn into transformations of moment matrices, or rather 
transformations of the operators that are associated with them. 

We first examine the classical existence problem for moments, culminating in 
a new proof of the existence of a Borel measure on R or C with a specified list 
of moments. Next, we consider moment problems associated with affine and non- 
affinc IFSs. Our main goal is to determine conditions under which an intertwining 
relation is satisfied by the moment matrix of an equilibrium measure of an IFS. 
Finally, using the famous Hilbert matrix as our prototypical example, we study 
boundedness and spectral properties of moment matrices viewed as Kato-Friedrichs 
operators on weighted spaces. 
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ABSTRACT 



Preface 

Moments of Borcl measures /i on M'^ have numerous uses both in analysis and 
in applications. For example, moments are used in the computation of orthogonal 
polynomials, in inverse spectral problems, in wavelets, in the analysis of fractals, 
in physics, and in probability theory. In this paper, we study some well-known and 
perhaps not-so-well-known aspects of moment theory which have been motivated 
by the study of both the classical literature on moments and the newer literature 
on iterated function systems (IFSs). 

Over the last hundred years, since the time of Lebesgue, mathematicians have 
adopted two approaches to measures on a topological (Hausdorff) space X. In 
the first, measures are treated as functions on some sigma-algebra of "measurable" 
subsets oi X . In the second, measures are realized as positive linear functionals 
on a suitable linear space C of functions on X. If we take C to be the compactly 
supported continuous functions of X, then Riesz's theorem states that the two 
versions are equivalent. Starting with a measure /i on the Borel sigma algebra, 
integration with respect to fi yields a functional L defined on C, and conversely (by 
Riesz), every positive linear functional L on C takes the form of integration against 
some Borel measure /i. The conclusion from Riesz's theorem asserts the existence 
of fl. 

The moment problem is an analog, taking X to be M.'^ and replacing C with 
the linear space of polynomials in d variables. The issue, as before, is to construct 
a measure from a functional. Since the polynomials are spanned by the monomials, 
a functional L is then prescribed by a sequence of moments. Hence, we have the 
moment problem: determine a measure /i from the sequence of moments. Given 
the moments, one can ask about the existence of a measure having those moments, 
the uniqueness of the measure if one exists, and the process for constructing such 
a measure and determining some of its properties. We will touch on each of these 
topics in this Memoir. 

In Chapter [1] we introduce our notation, definitions, and conventions. We have 
collected here for the reader's convenience some of the tools from operator theory 
and harmonic analysis which we will use throughout. In Chapter [2l we review the 
classical moment existence problem, which asks the following question: Given a 
list of moments {mk}k>o, does there exist a Borel measure /i such that the fc*'' 
moment of /i is rrifc? Alternately, we can arrange the moment sequence in a Hankel 
matrix in order to apply operator-theoretic techniques to this problem — a theme 
we carry throughout the Memoir. Using tools of Kolmogorov, Parthasarathy, and 
Schmidt, we provide a new approach to this old problem, which was first settled by 
Riesz (see iRie23. ) nearly 100 years ago. Our main result here is a new proof for 
the moment existence problem in both R and C. We return to this same theme — 
the classical moment uniqueness problem — in Chapter [51 There, we describe in 
detail the theory of self-adjoint extensions of symmetric operators, which helps us 
determine precisely when a list of moments has more than one associated measure. 

In Chapters [3] and IH we explore the difficult problem of computing moments 
directly for an equilibrium measure arising from an iterated function system (IFS). 
An IFS is a finite set of contractive transformations in a metric space. A theorem 
of Hutchinson [HutSlj tells us that for each IFS, there exists a unique normalized 
equilibrium measure n which is the solution to a fixed point problem for measures. 
We show that every IFS corresponds to a non-abelian system of operators and a 
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fixed point problem for infinite matrices. We then prove that the moment matrix 
M = M*^^-* is a solution to this matrix fixed point problem, and in turn we exploit 
this fixed point property to compute or approximate the moments for /i in a more 
general setting than the affine cases studied in |EST06j . 

As shown in |EST06j . it is not straightforward to compute the moments for 
even the simplest Cantor equilibrium measures, but we can compute the moments of 
an equilibrium measure /i directly in the affine IFS case. However, we generally have 
no choice but to approximate the moments in non-affine examples. In particular, 
our results can be applied to real and complex Julia sets, in which there is much 
current interest. The non-affine moment approximation problem is surprisingly 
subtle, and as a result, we turn to operator-theoretic methods. Affine IFSs are 
considered in Chapter [Sj while non-affine IFSs and operator theory are considered 
in Chapter [31 In addition, there are associated results about spectral properties of 
moment matrices for equilibrium measures in Chapter [71 

Infinite Hankel matrices cannot always be realized directly by operators in the 
l"^ sequence space. We have been able to apply operator theoretic results more 
widely by allowing matrices to be realized as operators on a renormalized Hilbert 
space when necessary. We turn to this problem in Chapter[5l where we introduce the 
operator theoretic extensions of quadratic forms by Kato and Friedrichs [KatSOj . 
The quadratic form we use in this context is induced by the moment matrix M^'^\ 
Using the Kato-Friedrichs theorem, we obtain a self-adjoint operator with dense 
domain. This generally unbounded Kato-Friedrichs operator can be used to obtain 
a spectral decomposition which helps us understand the properties of the quadratic 
form which gave rise to the operator. Often, renormalized or weighted spaces allow 
us to use the Kato-Friedrichs theorem in greater generality. 

We continue to use Kato-Friedrichs theory in the remaining chapters to under- 
stand the spectral properties of the moment matrix operator. In Chapter [6l we use 
the classical example of the Hilbert matrix and its generalizations from Widom's 
work [Wid66| in order to explore spectral properties of the moment matrix A/*^^) 
for general measures. In particular, we show that the moment matrix is unitarily 
equivalent to a certain integral operator. In Chapter[71 we further explore the spec- 
tral properties of the moment matrix and present some detailed examples. Finally, 
Chapter [51 uses spectral theory as a tool to reexamine the classical moment prob- 
lem, this time considering not only the existence of a measure having prescribed 
moments, but also the uniqueness. 

Readers not already familiar with the theory of moments may find some of 
the classical references useful. They treat both the theory and the applications 
of moments of measures, and they include such classics as jST43j . [Sho47| . and 
[Akh65| . The early uses of moments in mathematics were motivated to a large de- 
gree by applications to orthogonal polynomials [Sze75| . More recent applications of 
orthogonal polynomials are numerical analysis, and random matrix theory, random 
products, and dynamics. These applications are covered in [Lan87b| . [Lan87a] . 
[Dei99| . while the edited and delightful volume [Lan87a| includes additional ap- 
plications to geometry, to signal processing, to probability and to statistics. 
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CHAPTER 1 



Notation 

In this section, wc introduce some notation, conventions, and definitions needed 
throughout the paper. 

1.1. Hilbert space notation 

We use the convention that the inner product, denoted {u\v), is hnear in the 
second position v and conjugate Hnear in the first position u. This choice results 
in fewer matrix transpositions in the type of products we will be computing and 
therefore yields more straightforward computations. 

We will use Dirac's notation of "bras" and "kets" . For vectors u and v in a 
Hilbert space, the inner product is denoted "bra-ket" {u\v). In contrast, the "ket- 
bra" notation denotes the rank-one operator which sends a vector x in the 

Hilbert space to a scalar multiple of u. Specifically, it is the operator given by 
X I— > {v\x)u. The Dirac notation for this operation is 

\u) (v\ x) = {v\x)u, 

or also is sometimes written to preserve the order \u) {v \ x) = u{v\x) . 

If we consider these Hilbert space operations from the point of view of matrices 
and Euclidean vectors, we denote a column vector by which we call a "ket". 
Similarly, we denote a row vector by {u\ and call it "bra" . The notation for the inner 
product and rank-one operators above now are consistent with the actual matrix 
operations being performed. Further, the algebraic manipulations with operators 
and vectors are done by simply merging "bras" and "kets" in the position they 
naturally have as we write them. 

1.2. Unbounded operators 

We state here some definitions and properties related to unbounded operators 
on a Hilbert space. For more details, see [RS801 ICon90| . 

Definition 1.1. An operator F which maps a Hilbert space Hi to another 
Hilbert space TL2 is a linear map from a subspace of TLi (the domain of F) to 7^2- 

It will be assumed here that the domain of F is dense in Tii with respect to 
the norm arising from the inner product on Tii. If F is not a bounded operator, 
we call it an unbounded operator. 

Definition 1.2 f [Con90| ). An operator F from Hi to "^2 is closed if the 
graph of F, {(x. Fx) C Tii x H2 ■ x G dom(_F)} is a closed set in the Hilbert space 
Hi X H2 with inner product ((xi, X2)|(yi, 2/2))-Hi xWs {xi\yi)Hi + {x2\v2)n2- If 
there exists a closed extension to the operator F^ then F is called closable. In that 
case, there exists a smallest closed extension, which is called the closure of F and 
is denoted F. 
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Definition 1.3. Let F be an unbounded operator with dense domain dom {F) 
from Til to 7^2- Let D be the set of all y E H2 such that there exists an a: G Hi 
such that 

{Fz\y)n2 = {z\x)h, 

for all z E dom(i^). We define the operator F* on the domain D by F*y = x. 
Note that x is uniquely determined because F is densely defined. F* is called the 
adjoint of -F. 

It follows from these definitions (see [RS80| ) that F is closable if and only if 
F* is densely defined, and that F* is always closed. 

Definition 1.4. An operator F on a Hilbert space H is symmetric if dom(F) C 
dom(F*) and Fx ~ F*x for all x E dom(F). F is self-adjoint if F is symmetric 
and dom(F) = dom(F*); i.e. F* = F. 

Often, the term hermitian is also used for a symmetric operator. We see that a 
symmetric operator F must be closable, since the domain of F* contains the dense 
set dom(F), and is therefore dense. 

Definition 1.5. A symmetric operator F is essentially self-adjoint if its closure 
F is self- adjoint. 

If an operator F is bounded on its dense domain and symmetric, it is essentially 
self-adjoint. 

1.3. Multi-index notation 

In W^, for d > 1, we will need to use multi-index notation. Here, a and (3 
denote multi-indices belonging to Nq, where Nq is the Cartesian product 

(1.1) Ng^ Np X ■ ■ ■ X Mq . 

d times 

Following standard conventions, the sum a + /? is defined pointwise: 

(1.2) a + /3 = (a,+A)ti- 

Using this notation, we have the following integral expression: 

(1.3) ( x°'dfi{x)^ [ x^'x^^ ■■■x'^-'dti{x). 

1.4. Moments and moment matrices 

Let X be one of M, C, M'^ with Borel measure /x. When X = M, we define the 
^th Qj.(;[gj. moment with respect to /i to be 

(1.4) a;M/i(x). 

Jr 

If the moments of all orders are finite, we will denote by M^^-* an infinite matrix 
called the moment matrix. The moment matrix in the real case has entries 

(1.5) M^^^ = m,+, = J x'+^ dfi{x). 

Throughout this monograph, the real moment matrices will be indexed by Nq x Nq 
- in particular, the row and column indexing both start with 0. 
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Definition 1.6. An infinite real matrix M whose entries are indexed by No xNo 
is called a Hankel matrix if 

for all values of k for which these entries are defined. 

The moment matrix defined in Equation (|1.5p is a Hankel matrix. This justifies 
our notation above referencing the (i, j)"^ entry of M^^^ with i + j. The following 
examples will be used throughout the paper to illustrate our techniques and results. 

Example 1.7. Lcbesgue measure on [0, 1]. 

Let fi be the Lebesgue measure supported on [0, 1]. Then the moment matrix 
for is 



1 1 



This matrix is often called the Hilbert matrix. We will examine the properties of 
the Hilbert matrix in greater detail in Section 16.11 

Example 1.8. Dirac point mass measure ji = 5i. 

Let /i be the Dirac probability mass 5i on R. Then the moments for /i are 
Af/^^) ^ j a;»+J' d^ = = 1, 

so the moment matrix entries are all 1. In this case, we see that the moment matrix 
does not represent a bounded operator on l'^ . This example will be studied further 
in Chapters [5] and [71 

Example 1.9. The measure = e""^ dx. 

Let ^ be e~^dx on = (0,c»). Using integration by parts, a quick induction 
proof shows that the moments for ^ are 

Again, we see that the moments increase rapidly as i,i increase, so A/'^' cannot 
be a bounded operator on We will discover more properties for this matrix in 
Chapters E] and m 

Example 1.10. Measures with moments from the gamma function. 

Let 



T{k) = / s'^-^e-Ms = (fc - 1)! 



when k G Z+. Set ji = e ^ ^ c\x. is a Gaussian measure with support on R. We 
have the odd moments equal to zero, and the even moments are given by 



for k E Z+. 

If X = C, the moments are now indexed by No x No and are given by 
(1.6) my = / z'z^dfi(z). 
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If every moment is finite, then the complex moment matrix M^'^^ is given by 

(1.7) A/g^ =my ==£z'z^dAi(z). 

Notice that a complex moment matrix will not in general have the Hankel prop- 
erty that arises in real measures. We also observe that the complex moments are 
equivalent to the inner products of monomials in the Hilbert space i^(/i): 

(1.8) m,y = / z'z=d^l{z) = 

Jc 

Moments and moment matrices arise naturally in the study of orthogonal polyno- 
mials in i^(/i). 

Example 1.11. If /i is the Lcbesgue measure supported on the unit circle T in 
C, then the moment matrix is the identity matrix: 

m\^^ = [ z^z'' d/i = / e2'^*('=^J> dx = 5jk. 
Jt Jq 

As we noted above, the identity matrix is not Hankel. 



In the case where X = M.'^, we index the moments using the multi- index notation 
defined in Section [T31 Given a € Nq, we have the a-moment 

(1.9) [ x°'dn{x). 

The moment matrix for K'' is actually indexed by Nq x Nq and has entries 

(1.10) AM = mc,+f3= f x^+I'dfiix). 

When M has real entries indexed by Ng x Ng, we will also call M a Hankel 
matrix if 

for all 7 for which these entries are defined. We see from Equation (|1.10p that the 
moment matrix for a measure on M'' {d > 1) also has the Hankel property. 

1.5. Computations with infinite matrices 

In a number of applications throughout this paper, we will have occasion to use 
infinite matrices. They will be motivated by the familiar correspondence between 
linear transformations and matrices from linear algebra. Given a linear transfor- 
mation T between two Hilbert spaces, assumed to be infinite dimensional, then a 
choice of ONBs in the respective Hilbert spaces produces an infinite matrix which 
represents T. Moreover a number of facts from (finite-dimensional) linear algebra 
carry over: for example composition of two transformations (and a choice of ONBs) 
corresponds to multiplication of the associated two infinite matrices. Moreover by 
taking advantage of the orthogonality, one sees in Lemma 11.131 that the matrix 
multiplication is convergent. 

Given an infinite matrix M, it is a separate problem to determine when there 
exists a linear transformation T between two Hilbert spaces and a choice of ONBs 
such that M represents T. This problem is only known to have solutions in special 
cases. We will address this further in Chapters [5] and [H We will encounter infinite 
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matrices which cannot be reahzed directly by operators in the ^^-sequence spaces, 
but for which an operator representation may be found after a certain renormahza- 
tion is introduced. 

We use the foHowing standard language and notation: if G = Gi,j is an infinite 
matrix indexed by No x No, and x is an infinite vector indexed by No, we will define 
the matrix-vector product Gx componentwise using the usual rule, provided the 
sums all converge. 

{Gx)i = ^ GijXj. 

Similarly, for infinite matrices G and H, the matrix product GH is also defined 
componentwise using the same summation formula used for finite matrices, provided 
these infinite sums all converge in the appropriate sense. Specifically, the («, j)*'^ 
entry of the matrix GH is given by 

(1-11) mkj = E g^mh,^,. 

k<£No 

As an aside, there is a parallel notion for such products which applies to matri- 
ces indexed by Nq x Nq and vectors indexed by Ng. Given M indexed by Ng x Nq, 
a vector c, and given a € Ng, we have 

(1.12) (Mc)„ = E Af„,^c^, 

provided that this sum converges absolutely. We particularly need absolute con- 
vergence here so that the sum can be appropriately reordered to sum over Ng. 

Definition 1.12. When the formal summation rules for a product of infinite 
matrices or a matrix-vector product yield convergent sums for each entry, we say 
that the matrix operations are well defined. 

Using our matrices G and H above, if the sums X^fceNo ^i.k^kj converge for 
all i,j £ No, then the matrix product GH is well defined. 

If G and H are operators on a separable Hilbert space, we can determine nec- 
essary conditions on their corresponding matrix representations so that the matrix 
products are well defined. If we take the Hilbert space to be 7i = £^(No) and let 
{cj : j G No} be the standard orthonormal basis in £^(No), i.e., ej{k) = Sj^k for 
j, fc S No, then we must first assume that G and H are defined on a dense domain 
which includes {ej}. This allows us to write the matrix representations of G and 
H . From there, we have the following result. 

Lemma 1.13. Let G and H be linear operators densely defined on £^(No) such 
that Gej, G*ej, and Hcj are well defined and in £'^(No) for every element of the 
standard orthonormal basis {ej}j^jqo- Then G = (Gij) and H = (Hij), the infinite 
matrix representations of G and H respectively, are defined and the matrix product 
{{GH)ij) is well defined. 

Proof. Recall that the matrix representation of the operator G is Gij = 
{ei\Gej)2 for i,j g No, provided Gej is in i"^ so that the inner product is finite. Our 
hypotheses imply, then, that the matrix representations of G and H exist. 

Given an operator G, the adjoint operator G* satisfies {G*u\v) = {u\Gv) for all 
u, V which are in the appropriate domains and for which G*u and Gv are in £^(No). 
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We now check for absolute convergence of the matrix product sums from Equation 

dniD: 

oo OO 



fc=0 fe=0 
oo 



J2\iG*e^e,){e,\He,)\ 



k=0 

< KG* e.|e,)|M [Y.\{ei\He,)n 

= ||G*ej||2!|^^ej||2 < oo. 

We use the Cauchy-Schwarz inequahty above since we know the vectors G*ei and 
Hcj are in i'^. Note that if G, H are bounded operators, the sum above is bounded 
by l|G*l|op||-ff||op for any choice of i,j G Nq. 

Once we know the sums from Equation (jl.lip are absolutely convergent, they 
can be computed: 



^G,,fci?fc,, ^ Y.^e,\Gek){ek\Hej) 



k=0 k=0 



J2{G*e.,\ek){ek\He, 



k=0 

= {G*ei\Hej) by Parseval's Identity 
= {eMGH)e,) 
= {GH),,,. 

□ 



Note that the conditions in Lemma 11.131 imply that the operator G is auto- 
matically closable. It is an immediate corollary that bounded operators satisfy the 
hypotheses of Lemma I1.13[ and thus their matrices will always have well defined 
products. 

Corollary 1.14. If G and H are bounded operators on £^(No), then the prod- 
uct of their matrix representations is well defined. 

In some of the applications which follow - for example the cases for which both 
G and H are lower triangular matrices - the range of the summation index will be 
finite. (See, for example, the proof of Lemma [3TT1 ) But we will also have occasion 
to compute products GH for pairs of infinite matrices G and H where the range of 
the summation index is infinite. In those cases, we must check that the necessary 
sums are convergent. 

There are two approaches to working with the product GH of infinite matrices 
G and H. One is the computational approach we have described above, and the 
other involves associating the matrices to operators on appropriate Hilbert spaces 
(see Chapters [5] and [6]). For many applications, the first method is preferred. In 
fact, there are no known universal, or canonical, procedures for turning an infinite 
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matrix into an operator on a Hilbert space (see e.g., |Hal67| and |Jor06j ). so often 
the computational approach is the only one available. 

1.6. Inverses of infinite matrices 

Computations in later sections will require the notion of inverse for infinite 
matrices. Let G be an infinite matrix indexed, as discussed in Section 11.31 by the 
set Nq X Nq. We will now make precise our (admittedly abusive) use of the notation 
G^^ for an inverse. First, in order to discuss computations, the index set for rows 
and columns must be equipped with an order. If c? > 1, we will use the order 
of the set Nq along successive finite diagonals. Our goal is to find an algorithm 
for the entries in the infinite matrix we denote G^^ . This does turn out to be 
possible for the infinite matrices we will be using in our analysis of moments and 
of transformations. 

Definition 1.15. Let E be an infinite matrix. We say that E is an idempotent 
if the matrix product E'^ is well defined and if E^ = E. 

Definition 1.16. Let G,H,Ei and E2 be infinite matrices. Assume that the 
matrix products GH and HG are well defined. We say that G is a left inverse of H 
if there is an idempotent matrix Ei such that GH = Ei. G is called a right inverse 
of H if there is an idempotent E2 such that HG = £'2 • If G is both a left and right 
inverse of iJ, G is called an inverse of H. 

While this notion of inverses for infinite matrices is not symmetric, the following 
lemma does justify the use of the term "inverse" . 

Lemma 1.17. Let G and H be infinite matrices such that both matrix products 
GH and HG are well defined. Suppose there is an idempotent Ei such that GH = 
El and EiG = GEi = G. Then the infinite matrix HG is an idempotent, denoted 
E2, which satisfies E2H = HEi and E2E1 ~ E2- 

Proof. First, we see that HG is idempotent: 

HGHG = H{Ei)G = HG. 

The formulas are also readily verified: 

E2H = HGH = H El 

and 

E2E1 = HGEi = HG — E2. 

□ 



A result of Lemma 11.171 is that if we know the matrix product GH = Ei is 
an idempotent such that EiG = GEi = G, then by Definition II. 161 G and H are 
inverses. 



Example 1.18. The following matrices arise in Example 17.171 with respect to 
measures which are convex combinations of Dirac masses. 



We are given the matrices 
1 



G = 



-1 




and iJi = 



8 
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We see that GHi = E where E is the idempotent (m fact, projection): 

"10 



E = GHi 



1 




Therefore, Hi is an inverse of G. Note that this inverse is not unique. The matrix 

1 



Ho 



1 f 1 

2 2 2 

111 
2 2 2 







also satisfies GH2 = E hence is also an inverse of G. 







CHAPTER 2 



The moment problem 

In this chapter wc introduce the moment problem, in which we describe the 
various positivity conditions that an infinite matrix M must satisfy in order to 
imply that there is a Borel measure fi on an ambient space X such that M is the 
moment matrix of /i, i.e. M = M^^^. These conditions wiU be different in the cases 
where X is a subset of R, C, or R'' for d > 1. We also examine whether or not 
the measures found will be unique. The existence of the measure ^ is discussed in 
Section [2.11 ^^i^d a procedure by Parthasarathy using a Kolmogorov construction is 
used in Section [2.21 to discuss a more concrete construction of fi. The connections 
of the moment problem to other areas of interest to mathematicians, physicists, 
and engineers are mentioned in Section [2.41 



2.1. The moment problem M = M^''^ 

Given an infinite matrix M , the moment problem addresses whether there exist 
measures fj. such that M = M^^\ There are known conditions on M which ensure 
that such a measure fj, exists. Our presentation will be brief and we will omit 
the proofs of the results stated here, as they are available in the literature, albeit 
scattered within a variety of journal articles. Fuglede's paper |Fug83| offers a very 
readable survey of the literature on moments in several variables, up to 1983. The 
two papers [BD061 IBD07| emphasize the complex case and include new results. 

Our intention in this section is primarily to give the definitions of the positive 
semidefinite properties on the infinite matrices in the real and complex cases. We 
will only treat the existence part of the moment problem, but there are a variety of 
interesting uniqueness results in the literature as well. Here we shall only need the 
simplest version of uniqueness: these are the cases when the measures fi are known 
a priori to be compactly supported. In this case, uniqueness follows as a result of 
an application of the Stone- Weierstrass Theorem. An example of nonuniqueness is 
seen in Section [5751 

We begin with the real case M'', for d > 1. Let V be the space of all infi- 
nite sequences c which are finite linear combinations of elements of the standard 
orthonormal basis {ea}a<£N^- Clearly I? is a dense subspace of the Hilbert space 
£^(Nq), and matrix- vector products are well defined on vectors in V. 

One of the ways to test whether a given sequence, or a given infinite Hankel 
matrix, is composed of the moments of a measure fj, is to test for a positive semidef- 
inite condition. While there are several such conditions in the literature, we isolate 
condition ([2.1^ as it summarizes a variety of features that will be essential for our 
point of view. Note that (|2.111 entails a separate verification for every finite system 

9 
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of numbers, hence it amounts to checking that all the finite truncated square ma- 
trices have positive spectrum. This in turn can be done by checking determinants 
of finite submatrices. 

Definition 2.1. A matrix M with real entries indexed by Nq x Nq is said to 
be positive semidefinite if 

(2.1) {c\Mc)i, = 51 H c^M^^pCp > 

aGNg/36Ng 

for aU ceV. 

We say that a given infinite matrix M is a (Nq, M'*)-moment matrix if there 
exists a positive Borel measure n on E** such that 

(2.2) M^,f, = Mi^l = I x"+'^d/i(a;) 

for all a,f3 Ng. Note that, as we observed in Section FOl if M is a moment 
matrix, it must have the Hankel property. 

The solution to the existence part of the moment problem in the M'* case is given 
in a theorem due to M. Riesz |Rie23| for d = 1 and Haviland [Hav35l IHav36j for 
general d. 

Theorem 2.2 (Riesz, Haviland). For any positive semidefinite real Hankel 
matrix M = (M^^^) indexed by Nq x Nq, there is a positive Borel measure ji on 
having moment matrix M . 

If M is an infinite matrix with complex entries, we say M is a complex moment 
matrix if there exists a positive Borel measure on C such that 

(2.3) M,,, = A'^? = / z'z^'^ii{z) 

for all 2, j G Nq. We define a property on the complex matrix M which is stronger 
than the positive semidefinite property. 

Definition 2.3 ( |BD06l IBD07| ). We say a complex infinite matrix M in- 
dexed by Nq X Nq has Property PDC if given any doubly-indexed sequence c = {ci,j} 
having only finitely many nonzero entries, we have 

(2.4) ^ Ci^jMi+t,]+kCkj > 0. 

Note that if we consider the special case where j,£ ^ in Definition [2?3l we get 
exactly the positive semidefinite condition on M given in Definition 12.11 

Theorem 2.4 ( |Fug83| , IBD061 lBD07p . If M is an infinite complex matrix 
indexed by Nq x Nq which satisfies Property PD<C, then there exists a positive Borel 
measure ^ on C such that M = M^^) . 

Remark 2.5. The condition in Theorem 12.21 that the matrix entries in M are 
real numbers cannot be dropped. To see this, let d = 1 and take ^ be a fixed nonreal 
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number. Let AI = (^"'^'"')j_fcgNo • We see that 

j k 




> 

for all c G v. The infinite matrix M is therefore positive semidefinite. Without the 
condition that the components of M be real, then, we would conclude that there is 
a measure on R with Af as its moment matrix. We can also verify, however, that M 
has Property PDC from Theorem l2.4l Let c = {cij} be a doubly indexed sequence 
in V (i.e. c has only finitely many nonzero entries). Then 



> 




This leads to the conclusion that there is a measure on C having moment matrix 
M. In fact, given ^ G C, we note that the Dirac measure 5^ at ^ has the matrix M 
as its moment matrix, since the entries of M are the evaluation of the monomials 



z-' z 



^ at i: 



It is shown in |E\tg83| that measures with compact support are uniquely de- 
termined by their moment matrices. The Dirac measure, therefore, is the unique 
Borel measure on C having moment matrix M . If we have ^ S C \ (i?), then even 
though M satisfies the positive semidefinite condition from Definition 12.11 there 
cannot be a Borel measure /i on R such that M = M^^K 

In some sense, what we have shown is that there is a measure, but its support 
is not restricted to M. One of the classical moment questions is whether one can 
determine the support of a measure ^ from the geometric properties of its moment 
matrix M^^'\ 



The simplest moment situation arises when the measure ^ = i5j is a Dirac mass. 
We will revisit this example frequently in the remainder of the Memoir. Observe 
that the infinite moment matrix for has rank one. We will later show that if fi 
is a finite convex combination of Dirac masses, then the associated moment matrix 
will be of finite rank. 
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2.2. A Parthasarathy-Kolmogorov approach to the moment problem 

The traditional existence proofs of Theorems 12.21 and 12.41 are not constructive. 
In this section, we carefuUy explain how the Parthasarathy-Schmidt theorem, The- 
orem 12.71 below, can be applied to the moment problems from Theorems 12.21 and 
12.41 Theorem 12.101 and Corollary 12.111 are restatements of the moment problem 
in the real and complex cases, respectively. We use Theorem 12.71 to produce the 
measures which appear in Theorem 12 . 101 and Corollarv l2.11l In Corollarv l2.11l the 
condition PDC (Definition l2.3p appears in a more natural fashion. In addition, the 
proof of Theorem 12.101 clearly shows why the Hankel assumption for the infinite 
positive definite matrix M is essential in the real case. 

The use of the Kolmogorov ideas is motivated by our applications to iterated 
function systems (IFSs) which begin in Chaptcr[3l We will be interested in moments 
of IFS measures. A key tool in the analysis of these IFS measures will be infinite 
product spaces, precisely such as those that arise in the Parthasarathy-Schmidt 
construction. 

To simplify the main idea, we carry out the details only in the real case, and 
only for d = 1. The reader will be able to generalize to the remaining real cases in 
R'', d> 1. We conclude with the application to the complex moment problem. 

Definition 2.6. Let 5 be a set, and let M : S* x S* ^ C. We say that M is a 
positive semidefinite function if 

^^c7A/(s,i)ct > 

for all sequences {crises CE V. (Recall this means the sequences have only finitely 
many nonzero coordinates.) In the real moment problem, the matrix M is real, 
and we restrict to sequences {c^jsgs with entries in R. 

Theorem 2.7. ( [PS72| Theorem 1.2]) Suppose S is a set, and suppose M : 
S X S —f C is a positive semidefinite function. Then there is a Hilbert space Ti and 
a function X : S TL such that TL = sp{A(.s) : s G S*} and 

(2.5) M{s,t)^{X{s)\X{t))n. 
Moreover, the pair (7i, X) is unique up to unitary equivalence. 

Remark 2.8. We outline here two choices for the pair (7i, X) in the special case 
where S* = No and M is given by an infinite positive semidefinite matrix. The first 
pair is the the one constructed in |PS72j using Kolmogorov's extension principle 
on an infinite product space. This particular choice will allow us in Theorem 12. 101 
to express concretely the measure satisfying the moment problem in the special 
case where M is Hankel. 

The second pair is formed by constructing a Hilbert space from a quadratic 
form using the matrix M . We will be using similar completions later in the paper 
to work with moment matrices, so this is a natural approach. Note that Theorem 
12.71 tells us that our two Hilbert spaces are isometrically isomorphic. 

Proof 1. Let S be the natural numbers No and let Q, be the set of all functions 
from No into the one-point compactification M of R: 
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{fl must be compact in order to use a Stone- Weierstrass argument in the construc- 
tion.) 

We now describe the Gaussian construction and its Kolmogorov consistency. 
For now, assume that M is positive definite in the strict sense. Let J ~ {ii, . . . , ip} 
be a finite subset of Nq . The positive definite function M gives rise to a (strictly) 
positive definite matrix Mj which is formed by choosing elements from rows and 
columns in M indexed by J. Let Pj be the Gaussian measure on ilj := M with 
zero mean and covariance matrix Mj, such that Pj has Radon-Nikodym derivative 
fj with respect to Lebesgue measure, where fj is given by 



where we denote ujj = (uji-^ , . . ■ , ) . 

Note that Equation (|2.6p uses the strict positive definite condition because the 
inverse MJ^ is required. The density in Equation (|2.6p is the standard multivariate 
normal density for random variables X{ii), . . . , X(ip) with mean 0. 

By Kolmogorov's extension theorem |Kol50[ "Fundamental Theorem." p. 29], 
in order for the family of measures {Pj : J finite} to define a measure P on all of 
n, the measures Pj must be consistent. Suppose J and K are both finite subsets 
of No, where J d K . In this context, consistency means that if g is a function on 
f2j which is extended to a function G on Qk depending only on the variables in J, 
then 



If we consider this problem in the coordinates which diagonalize the matrix Mx, 
we see that the variables in K\J will integrate to 1. and the measures Pj are 
indeed consistent. Note that consistency in one set of coordinates does not imply 
consistency in another set of coordinates. 

If M is not strictly positive definite, then for any finite set J C No, we can 
change coordinates to diagonalize AIj and write flj = ker(Mj) © kcr(Afj)-'-. We 
take the measure Pj to be the Gaussian given by fj on ker(Mj)-'- and 6a, the 
distribution having mean and variance 0, on ker(Mj). Then the covariance matrix 
of Pj to is exactly the diagonalized Mj, and Kolmogorov consistency is maintained. 

By the consistency condition, the measures Pj defined on finite subspaces of 
can be extended to a measure P on fi. We now take the Hilbert space in Theorem 
12.71 to be 7i = i^(ri, P) and the map X : No ^ H to map n G No to the projection 
map onto the n**^ coordinate of the element a; e il: 



Then we see that H is the closed linear span of the maps {X{n) : n £ No} and 



We can now call M the covariance function. 

To show the uniqueness statement, suppose Hilbert spaces Hi and TI2 with 
corresponding functions Xi and X2 satisfy {Xi{s)\Xi{t))T-Ci = M{s,t) and Hi = 
sp{Xi(s) : s £ S} for i — 1,2. For each s € 5, let W be the linear operator such 
tha,tW{Xi{s)) — X-zis). defined this way is an operator, since if CiXi(si) = 



(2.6) 





X{n){Lu) = LLi{n). 
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then 



H2 



\ i i 

^c-Cj{X2{si)\X2{sj))n2 




= J2'CiCj{Xi{si)\Xi{sj))n, 



CiXi{st) 



= 0. 



■Hi 



By linearity and density, W extends to a map from Tii onto Ti.2, and by the asso- 
ciation of the inner products with the matrix M, we see that W is unitary. □ 

Proof 2. As in the first proof, let S = No- Let V be the set of all maps from S to 
C (i.e. complex sequences) such that only finitely many coordinates are nonzero. 
Given the positive semidefinite function Af , define the sesquilinear form 



S{vi,V2) 



sGNo tGNo 



This is clearly a sesquilinear form and therefore yields a quadratic form Q which is 
a seminorm 

Q{v) = \\vrM = S{v,v). 

The set NuIIm = {v £ : Q{v) = 0} is a subspace of T>, so the quotient space 
I?/Nu11m is an inner product space. We complete this space to form a Hilbert space 
which we denote Tig to emphasize the dependence on the quadratic form Q arising 
from the matrix M. We note here that this construction of a Hilbert space from 
a given positive semidefinite function M is unique up to unitary equivalence. We 
will make use of this Hilbert space completion of the quadratic form Q again in 
Chapters [5] and [H 

We next define the map X : No Hq by 

[xmt)-sAt)^[l l=\. 

Clearly (5s G I? C TCq for all s e Nq, and in fact, V is the linear span of {5s 
No}. Therefore, Hq = sp{X(s) : s G No}. We also see that 

{X{s)\X{t)) = '5MMM5t{v) 

u,v£No 

= M{s,t). 



s e 



□ 



Remark 2.9. There arc subtleties involved in these Hilbert space completions. 
In the pre-Hilbert space (before the completion), we will typically be working with 
a space T) of all finitely supported functions on S. In other words, T) is the space of 
all finite linear combinations from the orthonormal basis {5s}sgs for £"^{3). When 
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the completion to a Hilbert space is done, T) will be a dense linear subspace in 
the completion. Here, "dense" refers to the norm being used. In Chapter [5l for 
example, the norm is a weighted P norm. 

Caution: In these alternative Hilbert space completions, say Tig above, explicit 
representations of the vectors in TLq arc often not transparent. The useful geometric 
representations of limits of concretely given functions may be subtle and difficult. 
Notions of boundary constructions reside in these completions. In fact, this subtlety 
is quite typical when Hilbert space com pletion s are used in mathematical physics 
problems. Such examples are seen in 



JOOO, IJorOO] , where symmetries result 



in separate positive definite quadratic forms, and hence two very different Hilbert 
space completions. 

Wc arc now ready to state our existence result for the real moment problem, us- 
ing the language of Kolmogorov and Parthasarathy-Schmidt (the first proof above). 
Given a matrix M , our goal is to find a measure /i such that M is the moment ma- 
trix for /i. The previous theorem gives us a Hilbert space and a map (7i, X), unique 
up to unitary equivalence, which we can use under the right condition [M a Hankel 
matrix) to determine a solution to the moment problem. This solution will not, in 
general, be a unique solution. 

Theorem 2.10. Let M be an infinite matrix satisfying the positive semidefinite 
condition i2.1\) . and let Af(0,0) = 1. Let flnan be the measurable subset offl given 
by 

(2.7) ^Han =={ujen : Lu{k) = [uiO)f for all k e Nq}. 

Then Kolmogorov 's extension construction yields a probability measure Puan on 
^Han with M as its moment matrix if and only if M satisfies the Hankel property 
from Definition \1.6\ 



Proof. {=>) Suppose a measure Pnan exists as stated. Theorem 12.71 defines 
the map X such that X{k){uj) = uj{k) = [w(O)]'' for lo € ^Han- For the covariance 
function we have 



(2.8) M{j,k)^ / XU)X{k)dPHan- 

Using Equation ()2.7p . we then find 

M(j,k) = / [X{Q){u)V[Xmio)fdPHan{0j) 
X^x''diPHanOXiO)-^)ix) 

x^+''d^i{x) 



where we define the measure by 

(2.9) fi^ PHanoXiOyK 

Since this formula shows that M is a Hankel matrix, we have proved one implication. 

(<^) Conversely, if M is a Hankel matrix, then Kolmogorov consistency holds 
on the subset flHan C because i^Han is measurable. By the extension principle, 
we get a measure space {ftHan, Pnan) such that the measure fi = Pnan ° ^(0)""'^ 
on M is a solution to the moment problem for M — M^^^^ □ 
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There is an analogous result for the complex case, which we describe briefly. 
Let M : No X No ^ C be a function satisfying Property PDC, i.e., 

(2.10) CijM,+e,j+kCk,i > 

for all doubly- indexed sequences c = {cy} G V. Define an induced function M : 
X Ng ^ C by 

(2.11) Miii,j),ikJ))^M{t + l,j + k). 

One can readily verify that M satisfies the positive definite condition given in 
Theorem [O 

Corollary 2.11. When Theorem \2.7\ is applied to the function M above, we 
get a pair {X,TC) where the Hilhert space may he taken to be Ti. ^ L'^{<C,ij,), and the 
function can be given by X{i,j) = z^^^ for all {i,j) G Nq. 

Proof. Given these choices of H and X, we see that 

{X{t,j)\X{k,l))H^ / I^z^¥dM(z) 
Jc 

by Theorem [im Then, 

/ ;^z^+^d/i(z) = (z*+'|z-'"+'^)^.(^) 
Jc 

= M{i + l,j + k) =M{{i,j),{k,l)). 

This verifies our desired result and also shows that ^ is a solution to the complex 
moment problem for the matrix My = M[i,j). □ 

2.3. Examples 

Example 2.12. The measure ^ = e~^da;. 

We showed in Example II .91 that the measure e~^dx on the positive reals has 
moment matrix 

A<-' = (^ + j)! 

for all i,i G No. This measure does not have compact support, but it is an example 
of a case where the solution to the moment problem M = M^''^ is unique. We know 
this because the Laguerre system of orthogonal polynomials is dense in L^(/i,R+). 

Recall our earlier observation that this is also an example which illustrates 
that infinite Hankel matrices cannot always be realized directly by operators on 
£^-sequcncc spaces. However we will show that an operator representation may be 
found after a certain renormalization is introduced. This will be an example of a 
general operator theoretic framework to be introduced in Chapter [5] 

We next describe a series of examples of measures having moments involving 
the Catalan numbers. We will revisit these examples in Chapter [71 where we will 
be able to say more about the operator properties of the moment matrices. 

The k^^ Catalan number Ck is defined 
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where the first Catalan number, Co, is 1. We define Bk to be 

k 

The Catalan numbers satisfy the following relation: 

k 

Ck+1 = CnCk-n- 
n=0 

There is an explicit formula for the generating function associated with the Catalan 
numbers: 



(2.12) Gcat(x) = V Cx'^ = ^ 

2x 

k=Q 

The radius of convergence for Gcat (a-') is j ■ 

The generating function GBin {x) associated with the 's has the same radius 
of convergence; in fact, 

(2.13) Gcat(x) = - f GBiniy)dy, 

X Jo 

and 

(2.14) GBin(a:) 



VI -4a;' 



The Hankcl matrix M*-^**' = (Gj+fc) is positive definite because every principal 



^2 

submatrix has determinant 1. 



Example 2.13. Wigner's semicircle measure. 
The measure d^ is given on (—2,2) by 



A[i{x) = — Ax 

2n 

and is otherwise. A simple calculation shows that 

m2k = J x'^''dn{x) = Ck 

and 

m2k+i = J x^''+'dfi{x) = 0. 



Example 2.14. The secant measure. 
The measure d^i is given on (—2,2) by 



d/i(x) = — dx 
Try 4 — 



and is otherwise. In this case. 



J x^''d^i{x) = Bk 
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and 

f'2 



TO2fe+i = / a;2''+M/^(x) = 0. 



Example 2.15. The half-secant measure. 
The measure dfi is given on (0, 2) by 

2 

Try 4 — a;^ 



and is otherwise. In this example, both the even and odd moments are nonzero: 

m2k = / x'^''dfi{x) = Bk, 
Jo 



and 

7r(fc + l)V fc 



m2k+i = / .T2'^+M/^(a;) 



2.4. Historical notes 

Our positive semidefinite functions M in the form (|2.ip have a history in a 
variety of guises under the names positive semidefinite kernels, positive hermitian 
matrices, reproducing kernels, or positive definite functions, among others. In these 
settings, a positive semidefinite function is just a function in two variables, or, 
equivalently, a function M defined on S" x 5 where S" is a set. 

There is a broad literature covering many aspects of positive semidefinite ker- 
nels, especially in the case when 5' is a domain in a continuous manifold; S may 
even be a function space. Function spaces lead to the theory of reproducing kernels, 
which are also called Bergman-Shiffer-Aronszajn kernels. See jAro50j . [BS52| for 
a classical exposition and |BCR84] , |Dut04] , |DJ06| for a modern view. 

If S" is a group or a semigroup, and if M is a positive semidefinite kernel, we 
will adopt additive notation "-f" for the operation in S and restrict attention to 
the abelian case. Of special interest are the cases when the function M has one of 
two forms: 

(i) M{s,t) ~ Fi(s — t) for some function Fi on S 

(u) M{s,t) = F2{s + t), F2 a function on S. 
In the first case, we say that the function Fi is positive semidefinite on S. Equiva- 
lently, we say that M is a Toeplitz matrix. There is a substantial theory of positive 
semidefinite functions on groups (and semigroups), see e.g., |HR79j . |BCR84] : 
positive semidefinite functions play a central role in harmonic analysis. 

In the second case (ii), M is often called a Hankel matrix. In the following chap- 
ters, we develop operator theoretic duality techniques to study positive semidefinite 
kernels in the form of infinite Hankel matrices. 

(1) We define a (generaly unbounded) operator F which maps sequence spaces 
£2 or their weighted variants i'^{w) into function spaces L'^{fJ,). This crucial 
operator F : £'^{w) ^^(m) is defined on the dense subspace V of finite 
sequences in and sends a sequence into the generating function (or 
polynomial) . 
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(2) Given the Hilbert spaces £'^{w) and L^{^), we introduce a duality with 
the use of adjoint operators. If is a given operator from P{w) to 
its adjoint operator maps from L'^in) to £'^{w). 
As a byproduct of our operator analysis for the spaces £^{w) and L^ifJ.), we get 
operations on the two sides which are unitarily equivalent. The unitary equivalence 
is critical: it allows us to derive spectral properties of measures fi from matrix 
operations with infinite matrices, and vice versa. 

As a final historical note, we mention that because the operator F associates 
elements of sequence spaces with generating functions (or polynomials) , we initially 
obtained some of our results with formal power series in the spirit of Rota's um- 
bral calculus [RR78j , |Rom84] . We then found conditions (such as renormalizing 
sequence spaces) which guaranteed not just formal convergence but actual conver- 
gence. In Rota's umbral calculus, the "umbra" is a space of formal power series. 
In this Memoir, the "umbra" consists of the Hilbert spaces £'^{w) and i^(/i) and 
the closed linear operators F and F^ between them. For more about the book 
jRom84] and the umbral calculus's relation to other parts of mathematics, see the 
2007 review by Jorgensen [Jor07] . 



CHAPTER 3 



A transformation of moment matrices: the afRne 

case 

In this chapter we use the maps from an iterated function system to describe a 
corresponding transformation on moment matrices. In the case of an affine IFS, we 
show that the matrix transformation is a sum of triple products of infinite matrices. 

Given a measurable endomorphism r mapping a measure space {X, fi) to itself, 
we make use of the measure transformation 

(3.1) iit-^ flO T^^ . 

If {X, /i) and {X, jio t^^) have finite moments of all orders, we can also study the 
transformation of moment matrices 

(3.2) Af(M) ^ 

We are interested in the circumstances under which this transformation of 
moment matrices p.2|) can be expressed explicitly as a well-defined intertwining of 
the form 

(3.3) Af('^°^"') = A*M^>''^A 

for some suitable infinite matrix A. (We are using the notation A* here to represent 
the conjugate transpose of the infinite matrix A.) We then ask whether there are 
conditions under which A might be an operator on the Hilbert space (."^ . 

We begin by stating the result from [EST06| that the matrix A can be de- 
scribed explicitly in the case where r is an affine map of a single variable. We give 
a proof of this result in Section [3. II and then examine in Section [3. 21 the connections 
of this result to affine iterated function systems (IFSs) in one dimension. We also 
prove a stronger uniqueness result for the fixed point of iterations of the moment 
matrix transformation arising out of a Bernoulli IFS. We examine in Section 13.31 
the types of transformations that preserve infinite Hankel matrices, as a precursor 
to Chapter [4] in which we explore the existence of matrices A for more general 
measurable maps r. 

3.1. Affine maps 

We consider fjL a probability measure on M or C. We will consider the case 
where t is an afRne map, noting that a finite set {rf,} of affine maps on M or C can 
comprise an affine iterated function system (IFS). Such IFSs are used in the study 
of infinite convolution problems; see e.g., [Erd39l IJKS07bl IJKS08| . We begin 
with some definitions. 

Recall that the moment matrix M^'*) of /x is the infinite matrix defined by 
= /jj x''^^d^{x) when x is a real variable, and by mI'j = z^z^dn{z) when z 
is a complex variable. We naturally must assume that moments of all orders exist. 



21 



22 



3. A TRANSFORMATION OF MOMENT MATRICES: THE AFFINE CASE 



Certainly moments of all orders exist when ^ has compact support, but we will not 
always be restricted to compactly supported measures. Recall that the row and 
column indexing of M^'^' start at row and column 0. 

In the real case, we recall that every moment matrix il/''^' is a Hankel matrix. 
Our moment matrix transformation corresponding to r in this real setting must 
preserve the Hankel property, then, since M('"°^ is also a moment matrix. 

The following result is stated in [EST06] . Because one of our goals is to 
generalize the lemma, we include a proof here for completeness. This result gives a 
matrix A corresponding to an affine map r on C, so that the moment transformation 
is an intertwining by matrix multiplication. Note that the same matrix A works in 
the real case if we take t to be an afhne map of a single real variable. 

Lemma 3.1. jEST06[ Proposition 1.1, p. 80] Suppose t : C ^ C is an affine 
map (not necessarily contractive) given by 



t{z) ~ cz + b, 



where c,b G C. Let M^^^°'^ ^ be the moment matrix associated with the measure 
H o . We have 

(3.4) = A*M^^''>A, 

where A — (a^j ) is the upper triangular matrix with 



(i^c'fr'-^ i<j 
otherwise 



(3.5) Oi^j 

Proof. The {i,jY^ entry of the moment matrix A/'^""^ ^ is given by 



(3.6) 



= / {cz + bY{cz + by d^l{z). 



Even though A and M^'^^ are infinite matrices, the sum defining the («, j)*'' entry 
of the triple product A*M'^^'>A is finite because A is upper triangular. With this 
observation, the matrix product is well defined and we can compute the (i,J)*'' 
entry of the product A* Al'-'^'^A by 
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Taking advantage of the linearity of the integral, we compute the sums in k and £ 
to obtain 

= / (cz + by{cz + by dfi{z). 
Jc 

This gives us our desired result. □ 

3.2. IFSs and fixed points of the Hutchinson operator 

An iterated function system has an associated compact set X and a measure 
H supported on X, where both X and /i arise as unique solution to a fixed-point 
problem described in the following paragraph. A central theme in this Memoir 
is that every IFS in the classical sense corresponds to a non-abelian system of 
operators and a version of Equations (|3.7|1 and ()3.8|1 in which the moment matrix 
M = M^'^) will be a solution to an associated fixed-point problem for moment 
matrices. (See Proposition 13.31 and Sections 14. II and 15.41 ') One of the consequences 
is that we will be able to use the formula for M in a recursive computation of the 
moments, which is not easy to do directly for even the simplest Cantor measures. 

Let / be a finite index set. An iterated function system (IFS) is a finite collection 
{TijiG/ of contractive transformations on K and a set of probabilities {pi}i^i. Using 
Banach's fixed point theorem and the Hausdorff metric topology, it is proved in 
Hutchinson's paper |Hut81j that there is a unique compact subset X of M'', called 
the attractor of the IFS, which satisfies the equation 

(3.7) X^[jn{X). 

There is also a unique measure fj, supported on X which arises from Banach's 
theorem. 

Theorem 3.2 (Hutchinson, |Hut81j ). Given a contractive IFS {nj^fzi in R'^ 
and probabilities {pi}, there is a unique Borel probability measure p, on satisfying 

(3.8) M-Ep'(^°^r')- 

The measure /i is called an equilibrium measure. Moreover if pi > for all i Cz I , 
the support of p is the unique compact attractor X for {rijig/. 

First, let us examine an IFS on the real line M consisting of contractive maps. 
We will find that the moment matrix for the equilibrium measure of the IFS also 
satisfies an invariancc property corresponding to Equations (j3.7p and (|3.8p . 

Let / be a finite index set. Let a contractive IFS on M be given by the maps 
{Ti}i£i and the nonzero probability weights {pi}iei. Suppose fi is the invariant 
Hutchinson measure associated with this IFS. Then /i is a probability measure 
satisfying p.8p which is supported on the attractor set X. 

We now observe, as also noted in [EST06| . that the moment matrix for the 
equilibrium measure fi of an IFS on the real line also satisfies an invariancc property 
7?.(M(''') = A/(^) under the transformation 

(3.9) TZ : ^ M(E.e/P»('^°^r')). 
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In other words, M^''^ is a fixed point of the transformation TZ. We now state a 
uniqueness result. 

Proposition 3.3. Given a contractive IFS {Ti\iei on M and probability weights 
{Pi}iei with Pi € (0, 1) and^^^j Pi = 1, there exists a unique moment matrix 
for a probability measure fi which satisfies 

(3.10) 7^(M(^)) A/(^). 

This unique solution is exactly the moment matrix for the Hutchinson equilibrium 
measure for the IFS. Furthermore, for any Borel regular probability measure v sup- 
ported on the attractor set X of the IFS, TZ^{M^'^^) converges componentwise to 

Proof. We begin by proving the convergence of the moments. Let S be the 
transformation of measures supported on X given by 

(3.11) SM=J2p^i'^°r-'). 

From |Hut81| . we know that S{iy) maps probabiUty measures to probabihty mea- 
sures, and that for any Borel regular probability measure v supported on X, 
S^v — > ^ as n oo, where convergence is in the metric p which Hutchinson 
calls the L-metric on the space of measures |Hut81| . 

To define the L-metric, first recall the space of Lipschitz functions Lip(X, M) 
on a metric space X with metric d. Lip(X, M) consists of all functions (p : X ^M. 
such that there exists a constant C = < oo with 

(3.12) |(/.(x) - 0(j/)| < C4,d{x,y) for all x,y X. 
We then define the constant 

(3.13) i(0) = inf{C0 I C0 satisfies (j3T2|) l. 
and hence the space of functions 

Lipi = {0 e Lip(X,R) I L{(t>) < 1}. 
We can now define the L-metric: 

(3.14) p{p,,iy) = sup 1^ J <l)dfi~ J (Tidi^ : (/) e Lipi(X,R)| . 

When the measures are defined on compact spaces, the L-metric topology is equiv- 
alent to a weak topology, and therefore we have 

/ /d(5V)^ / fdp 
Jx J X 

for all functions / : X M which are bounded on bounded sets |Hut81| . 

Let /(x) = x^^^ , which is bounded on bounded sets for all choices of i, j G Nq. 
Inserting / into the convergence above yields the componentwise convergence of 
the moment matrix for S"'v to the moment matrix for p. 

Given that p satisfies the invariance property 

S{p) = /i, 

we next wish to conclude that 7?.(A/(^)) = Af(^\ We defined the transformation TZ 

by 
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Since Sji = /i, it is clear that 5/x and /x have the same moment matrices. Therefore 
7^(M(^)) = M^'^l 

It remains to be shown that the matrix M^''-' is the unique solution among 
moment matrices for probability measures on X to the equation TZ{M) = M . We 
have from |Hut81| that the measure /j, is the unique probability measure satisfying 

= fi. Suppose there exists another matrix M^''^ which is a moment matrix 
for a probability measure and which is a fixed point for TZ. Since 
Af('5''), the measures v and Sv have the same moments of all orders. By Stone- 
Wcicrstrass, since the measures have compact support and agree on the moments 
(hence the polynomials), the measures must be the same. Since Sv = v, we must 
have u ~ ji. □ 

Corollary 3.4. Let {Ti\iei be an ajfine contractive IFS on M, such that 
Ti(x) = CiX + bi. Let be probability weights and let Ai be the matrix given 

in Lemma \3.1\ which encodes ti for each i Cz I . The moment matrix transformation 

(3.15) M^'') ^ ^p.^U/t")/!, 

has a unique fixed point among moment matrices. Moreover, the fixed point is the 
moment matrix M^^^ for the Hutchinson equilibrium measure of the IFS. 

Given an affine IFS and its associated moment matrix transformation, we can 
state a uniqueness result for finite matrices which corresponds to and extends Corol- 
larv 13.41 Given any infinite matrix M, denote by Af„ the (n + 1) x (n + 1) matrix 
which is the upper left truncation [My]ij=o,i....,n of M. Let {t^I^Lq be a contrac- 
tive afhnc IFS on M, and let fi be the unique Hutchinson measure corresponding to 
the IFS and probability weights {pijjLg. Note that we have Ti{x) = CiX + bi, where 
Ci < 1 for each i = 0,l,...,fcin the most general affine IFS. 

We know from CoroUarv 13.41 that M^^^ is the unique infinite matrix which is 
fixed under the transformation in Equation (|3.15p . We now give the result that the 
truncated form of the transformation in Equation (|3.15|) has the truncated moment 
matrix as a fixed point, and moreover, that it is a unique fixed point among Hankel 
matrices. We remark here that some of the computations required for this proof 
are done in |EST06| . but they don't use them to state this uniqueness result. 

We begin by observing that, because the matrices Ai which encode afhnc maps 
are triangular (recall Lemma 13. ip . a truncation of the matrix product A*MAi is 
exactly equal to the product of the truncated matrices {Ai)'^Mn{Ai)n. 

Lemma 3.5. Given infinite matrices A,M such that A is upper triangular, 

(3.16) {A* MA),, = A^MnAr,. 

Proof. This follows immediately by the properties of block matrices assuming 
that the associated block products are well defined. In the matrices below, B and 
D are upper triangular: 



'B* 


" 




'E 


F' 




B 


c' 




'B*EB * 


C* 


D* 




G 


H 







D 







□ 

Given the affine IFS {Tij.f^g, {pij^^Q mentioned above, let Ai be the matrix 
given by Lemma l3. II which encodes each affine map Ti{x) — CiX + bi. Let ^ be the 
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unique Hutchinson equilibrium measure for this IFS. Using our truncation notation, 
define the matrix transformation TZn on (n + 1) x (n + 1) matrices M by 



(3.17) 



7^„(A/) =^p,(A*)„M(AO„. 



i=0 



We then have the foUowing result. 



Proposition 3.6. For each n e N, M ^ Mn^'^ is the unique Hankel matrix 
having Mq.o ~ 1 which is a fixed point for TZn- 



and Lemma [33] that for each n, mI''-' 



Proof. We know from Corollary 
a fixed point of TZn- We need only to prove the uniqueness. 

For n = 1, assume M is a 2 x 2 Hankel matrix with Mo,o 
7^l(M) = M. M is of the form 



IS 



1 such that 



M 



1 X 
X y 



Then, 



7^l(M) 



fe 

i=0 

fe 

i=0 





Ci 

1 



1 X 

X y 



bi + CiX 
1 



bi + CiX 
f- 2biCiX + 



xJ^PiCi + Yl,Pibi 

T.Pibi yJ2 Pi cf + 2xJ2 Pi bi Ct+J2 Pi 



1 X 

X y 



M. 



There is a unique solution for x and y in the equations x — x^piCi + "^Pibi 
and y = yY.Ptcf + 2xY.pib,Ci + Y.P^H-^ thus M = m[^'\ 

Assume next that is the unique n x n Hankel matrix which is a fixed 

point of TZn-i- Take M to be an (n + 1) x (n + 1) Hankel matrix with Mo,o = 1 
such that TZn{M) ~ M . By Lemma [23] and our hypothesis, M is of the form 



M = 



n-l 
gtr 



where B is n x 1. Due to the Hankel structure of M, all but one of the entries in 
B are known from mI^'' : 



B 



X 



It is straightforward to verify that, as in the n = 1 case, the terms in the 
matrix equation S{M) = M yield two linear equations in x and y which have 



unique solutions. Therefore M = M, 



(m) 



□ 
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The computations in Proposition 13. 61 to solve for x and y give a recursive equa- 
tion to compute the moments of an equihbrium measure in terms of the previous 
moments. 

Corollary 3.7 ( |EST06] ). Given the IPS as above with equilibrium measure 
1 / \ / \ 

= -— ^ — , E (7) [1) K^-^-'c^'M^ 

3.3. Preserving Hankel matrix structure 

When we work with affine IFSs for a real variable x, the moment matrix of a 
probability measure is a Hankel matrix whose (0,0)'^ entry is 1. We will examine 
here the sorts of matrix transformations of the form 

M A* MA 

which preserve Hankel structure, so that in the next section we can look for moment 
matrix transformations corresponding to more general maps t. Let Ti'^' denote the 
set of positive definite Hankel matrices whose (0, 0)*'^ entry is 1. 

We have already shown that the transformation TZ corresponding to an affine 
IFS given by 

(3.18) n{M) = E^''(A*^^^0 = Y.Pb^^^^°^''^ 

iei beB 

maps a moment matrix Af^^^ to another moment matrix Af^^'e^ P'^""^; ) . and there- 
fore maps Ti*^^) into itself. 

Definition 3.8. We say that the matrix A preserves Ti^^^ iiA*MA G 7^(1) for 
all matrices M en'^^\ 

In the following lemma, we refer to inverses of infinite matrices and to invertible 
matrices. A careful definition of these concepts is given in Definition II. 161 

Lemma 3.9. The following matrices preserve 

7 7^ 73 ..." 

1 27 372 ••• 

1 37 ••• 

1 



The inverse of D{S) is D{5 ^) (5 ^Q) and the inverse of G{'y) is G'(— 7). 

We are pleased to thank Christopher French for the proof of the following 
proposition. See also [FreOTj . which implicity uses Proposition 13.101 throughout, 
and the papers |SS06| and |Lay01| . 

Proposition 3.10. Suppose A = (cij) is an infinite upper triangular, invert- 
ible matrix which preserves Ti^^\ Then either A or —A is the product of matrices 
of the form D{S) and G{"f), where D and G are defined in Lemma \3.9l 



(3.19) D{S) = 



1 














6 














6^ 














6^ 



and G(7) 
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Proof. If A*MA e H^^^ for all M e Ti.'^^^ 
ao,o = — 1, replace A with —A. 
Denote A by 

"1 



then 



'0,0 



1, so ao = ±1- If 



A = 



ai,i ai^2 




0.2,2 



Since ai_i 0, multiply A on the right by D{l/ai^i). The upper 2x2 principal 
submatrix of AD{S) is 

1 ao,i/ai,i 
1 

Now set 7 = — ao,i/ai,i; the matrix AD{S)G{'^) now has the 2x2 identity matrix 
as its upper left principal submatrix. 

Now set A ~ AD{l/ai_i)G{—ao^i/ai,i); we know that A preserves Ti^^^ by 
Lemma [2121 We will show that A is the infinite identity matrix. We start with the 
k = 2 case (instead of A: = 1) to give more intuition. 

Let A be denoted 



(3.20) 



A 



1 ao 
1 ai 
a2 



Consider the upper left 3x3 principal submatrix of the matrix A* MA: 
(3.21) 



ao 



1 

mi 



mi ao + aimi+a2m2 
m2 * 
a2m2 * * 



bmce A* MA belongs to H'-^^ for aU M G n'^^\ we know that the above matrix is 
Hankel. We can treat mi and m2 as independent variables, so ao = oi = and 
a2 = 1. Therefore A actually has an upper 3x3 principal submatrix which is the 
identity. 

Continuing inductively, suppose the upper k x k principal submatrix of A is 
the k X k identity matrix and A has the form 



(3.22) 



10 

1 









ao 
ai 



1 afc_i * 
afe * 
* 



Multiplying A* MA, we find two expressions for the (fc, 0) entry: 
(3.23) mk = '^a.imi. 



i=0 
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Since this equation holds for aU choices of mi, . . . m^, we must have that ao = • • ■ = 

flfc-i = and ak = 1- 

Therefore, A = AZ?(l/ai^i)G(— ao,i/ai,i) = /, so we can write 

(3.24) A^Giao,i/ai,i)D{ai,i). 

□ 



CHAPTER 4 



Moment matrix transformation: measurable maps 

In the previous chapter, we focused on the case of afhne transformations com- 
prising IFSs, but there is a great deal of interest in concrete apphcations to nonafFine 
examples; e.g. real and complex Julia sets. We turn our attention to nonaffine mea- 
surable maps in this chapter. We begin with an arbitrary measurable map r on a 
measure space (X, /i), where X is a subset of M'^ for d > 1 or C, and the moments 
of all orders with respect to /i and /i o r^^ are finite. We ask whether the transfor- 
mation of moment matrices Af'^^ i— > M^f^°'^ ^ can be expressed as a matrix triple 
product A/t^""^ ) — A*M'^'^'>A, for A an infinite matrix. In the previous section, 
we stated the appropriate A when r is an affine map. We now seek to find an 
intertwining matrix A for more general maps r. 

We find that A, when it can be written down, need not be a triangular matrix. 
We therefore also will need to examine the hypotheses under which the formal 
matrix products we write down are actually well defined. 

4.1. Encoding matrix A for r 

The following analysis will be restricted to the real one-dimensional case. We 
let X C M be a Borel subset and let /i be a Borel measure with moments of all 
orders on X. Let r be a measurable map on X, so that ^ot~^ also is a measure on 
X with moments of all orders. We seek to find an infinite matrix A which enacts 
the moment matrix transformation from M*^^) to M^^""^ more precisely, such 
that 

(4.1) Af('^°^"') =: A*M^'''>A. 

In the space L^{fJ.), let V be the closed linear span of the monomials. (The 
monomials are functions since all moments are finite for /i.) The following results 
hinge on a careful description of the Gram-Schmidt process on the monomials, since 
the monomials could possibly have linear dependence relations among them in L^{fJ.) 
(See Example 14.11 below). Let {wjIjgNo be the set of monomial functions, but if 
is dependent on l,a;,a;^, . . . ,x^~^, we remove x'^ from the collection. In other 
words, Vj might not be , but it is a monomial function x^ for some k > j, and 
there are no finite dependence relations among the set {wjjjgN,-,. 

Example 4.1. Bernoulli IFSs which yield measures with linearly dependent 
monomials. 

Consider the afhne IFS on M of the form 



To{x) = Ax, Ti{x) — X{x + 2). 
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The properties of the Hutchinson equihbrium measure depend on the choice of the 
parameter A. If A < ^, the measure is supported on a generahzed Cantor set, while 
when A > ^ the measure is supported on an interval in the real line. 

We can transform these into measures supported on subsets of the unit circle 
T C C via the map x i— > e^^^". By [JP98| . when A = i, there is an orthonormal 
basis of monomials 

{z^ : fc = 0,1, 4, 5, 16, 17, 20, 21,...} 

for the corresponding measure contained in T. The collection of all the monomials, 
however, does not have finite linear dependence relations, as described in Remark 
14.21 On the other hand, if A = i, the corresponding circle measure does have finite 
linear dependence relations among the monomials. These cases also arise when 
A > i |JKS08| . 

Remark 4.2. A central theme in our study of moments in Chapters [3] and 
m is how the study of moments relates to the self-similarity which characterizes 
equilibrium measures for IFSs. Because we use operator-theoretic methods to study 
moment matrices, we encounter spectral information along the way. Historically 
(see }Akh65j ) the approach to spectral theory of moments in R went as follows: 

(1) Start with the monomials {x''}keNo viewed as a dense subset of (a subspace 

(2) Apply Gram-Schmidt to obtain the associated orthonormal polynomial 
basis {pfejfceNo- 

(3) The tri-diagonal matrix J representing multiplication by x in the ONB 
{pk}keNo (see Chapter [8]) gives rise to spectral information. 

However, the classical approach does not take into account the self-similar proper- 
ties that /i may have. Following [JP98| and |Jor06| . we can seek instead to encode 
the IFS structure directly into the analysis of moments. A case in point is the 
Cantor system mentioned above with scaling by j and two subdivisions (a measure 
fj. with scaling dimension i). It was shown in |JP98j that it is better to realize /i 
as a complex measure which is supported on the circle in the complex plane. The 
monomials we are then led to study are the complex monomials {z'^lfegN^ in L^(/i). 

This set of all the monomials has no finite linear dependence relations. To see 
this, we assume a polynomial p{z) = X]j=o '^j^'' supported on the circle is zero in 
L^{Ijl). Then for almost every point z in the support of ^, the polynomial is zero. 
Since this measure has no atoms (see |D J07j . [JKSOTa] ) . the complement of a set 
of measure zero must be infinite. We conclude that since p has infinitely many 
zeros, it is the identically zero polynomial and hence oq = oi = • • • = a„ = 0. 

The standard application of Gram-Schmidt on the monomials would produce 
an ONB of polynomials, but it would be different from an ONB of monomials 
such as {z^\k = 0, 1, 4, 5, 16, 17, 20, 21, ...}. Moreover, this approach misses the IFS 
scaling property of /i. 

Returning to our quest to encode the map r via a matrix we perform the 
Gram-Schmidt process on the finitely linearly independent monomials {"WjljeNo 
to construct an orthonormal basis {pfelfegNg of polynomials for V. Part of the 
definition of the Gram-Schmidt process gives that each polynomial pk is in the span 
of {vq, vi,. .. Wfc}, and hence is orthogonal to sp{vo, . . . , Wfc_i} = sp{po, • • -Pk-i]- 
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We can write the lower triangular matrix G which enacts Gram-Schmidt as follows: 

fc 

(4.2) ^Gfc,.«, -Pfe. 

Also, assume that is a probability measure, so po is the constant function 1, i.e. 
Po[x) = 1. 

In the cases where some of the monomials have been left out of the sequence 
{wj}, we define a moment matrix N^^^^ for /i to be 

(4.3) iVj,'^^' = {v,\vk)m^). 

This adjusted moment matrix A^(^) will be symmetric in the real cases but will not 
have the Hankel property. The moments contained in A^('^) are total in V. 

We place the condition on the map t that the powers are in the space V for 
all J G Nq. From here, we define the following transformations on V: 



(4.4) Rpk = vk 

(4.5) Tpk = Vk o T. 

R and T are well-defined operators on V since they are defined on an orthonor- 
mal basis. They might be unbounded, but their domains do contain the ONB 
elements {pfcjfceNo by definition. Therefore, we can express R and T in matrix 
form with respect to {pk}keNo- With a slight abuse of notation, we will also refer 
to these matrices as R and T respectively: 



(4-6) Rj.k = {P]\Rpk)mf,) = {pj\vk)mf,) 

(4-7) Tj^k = {Pj\Tpk)L^(j,) ^ {Pj\vkOT)L2(f,y 

Observe that the matrix for R is upper triangular, since each pj is orthogonal to 
sp{vo, . . .Vj-i} = sp{po, . . -Pj-i}. 

Lemma 4.3. The matrix R is invertible in the sense of DeHnition l 1 . 1^ and the 
matrix of R^^ is exactly the transpose of the matrix which enacts the Gram-Schmidt 
process on the monomials, i.e. 



Proof. We compute the matrix product RG^^ . Note that the triangular struc- 
ture of G*'' gives a finite sum, so the product is well-defined: 
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oo 



oo 



{pj\^G.j^kVk 



) 



k=0 




We need to show that this composition has a dense domain. One can readily verify 
that when a monomial Vj is expressed as a column vector with respect to {p^lfegNoj 
that column vector only has nonzero entries {pk\vj) for k < j. Moreover, by the 
upper triangular structure of G*'', the matrix- vector product G^^Vj must be in T). 
These are in the domain of the matrix for _R, so we can now conclude that the 
monomials Vk are all in the domain of the operator RG*^. Therefore, _RG*'' has 
dense domain and on that domain, i?G*'' is the identity. By Lemma [1.171 G**" and 
R are inverses of each other, and we write G''' = R~^. □ 

Next, we will discuss the adjoints R* and T* . These can be written down as 
matrices, but it is not always true that their domains are dense in V. We will show 
that there exists a renormalization of (/x) such that both R* and T* have dense 
domains in V. This is equivalent (see |Con90| Proposition 1.6, Chapter 10]) to 
saying R and T are closable operators. 

Given a set of nonzero weights w = {wijigNoi define the space Vw to be the set 
of all measurable functions which are in the span of the monomials with respect to 
the weighted norm 



To be precise, we see that / fpidfi must be finite for all i S Nq, but depending on 
the weights, observe that Vw could include functions which are not in L^{fJ,). The 
inner product, then, on Vw is given by 



fe 

If we consider the operator R as defined above, but as a map from Viu to V, 
then the adjoint i?^ : V Vw is given by 



on every g in the domain of (i.e. every g such that 7?* 5 e Vw) Observe that 
the weights change the adjoint, so we denote it by i?* . To verify Equation (|4.8p . 
we compute for / G Vw H Dom(i?) and g GV\ 



00 




(4.8) 




00 



k=0 
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which gives 



k=0 k=0 

We then compute 

^ (^fc l^^^'(^)^^ k)p,„ = II ( II ~ ^""^ \9)LH^.)Pk |PJ (Pj I/)l^(a.) 



The equahty of these two expressions verifies our definition of i?^, at least for the 
dense set of finite linear combinations of the orthogonal polynomials, hence holds 
for all / in the domain of R. 

In the next lemma, we produce a weighted norm on the space V using weights 
{wk} which guarantee that R will be closablc in the weighted space Vw 

Lemma 4.4. Given the measure space {X,fi) as above having finite moments of 
all orders and given the operator R defined in Equation li4-4^ , there exist weights 
{wijigNo such that R : Vw V is closable. 

Proof. Observe that Rl,g e Vw precisely when 
(4.9) ^«^fc|(Pfc|i?;5)L^(M)l' < oo- 

k 

Substituting in the definition of R*^ shows that Equation (|4.9p is true if and only if 



(4.10) 



< oo. 



Since /i has finite 0*'' moment, i.e. is a finite measure, we have L°°{^) C L^(/i). 
Set M.k J \ vk\dfJ.{x), and suppose / S L°°{pi). Thenf 



(4.11) 



(«fc|/)L2(,,; 



< 



WfW < ||/||L=o(^,)A^fc- 



Therefore, if {Ml/wk} € (\ then 

(4.12) L°"{^l) Cdom{R*). 

Wc thus explicitly define a choice of weights w = {wk}keNa by 

(4.13) Wk := {1 + k^)Ml. 

Since L°°{^) is dense in L'^{n), we now know that dom(i?J^) is dense in V, where 
w = {wk}keNo is defined in (|4.13p . Therefore R : Vw V is closable. □ 



The space Vw has the weighted orthogonal polynomials \ -76= > as an or- 

thonormal basis. We will denote these vectors {PfcjfeGNo- We next observe the 
following property of the matrix R!^R written in terms of the orthonormal basis 
{PfelfcsNo: 
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(4.14) KR^^, = mKRpDv^ 

= {Rpr\RpJ)LH,} 

{vi\vj)m,,) 



- *J 



So is a self-adjoint operator and is exactly a weighted version of the 

moment matrix. 

Next, we recall our definition (Equation (|4.5p ) for the operator T which maps 
Pi to Vi o T for each i € Nq- By our hypothesis on r, we know T is an operator on 
■p. We wish T to also be closable, so that T* is densely defined. As we discovered 
for R, this may require weights. Since the adjoint depends on the weights, we will 
denote it T^. 

A function / is in the domain of with respect to weights w = {uifcjfcgNo if 
T^f is in Vw The same computations used for replacing each Vk with Vk o r, 
show 



oo ^ 

We also find that g G dom(T^) if and only if 

oo ^ 

(4.16) V — (wfcOT|g)L2() 



Wk 



2 

< OO. 



The condition on r under which there are weights such that is densely 
defined on Vw is given below. 

Lemma 4.5. Let M'^ — \vk o T\d^i = / \vf;]d{fi o r~^). If M'^. < oo for all 
k € No, we define the weights w — {wfcjfegNo 

for which is densely defined from V to Vw 

Proof. Repeating the computation in Lemma wc find that / e L°°{ii) is 
in the domain of if and only if 



Y.—{M'nfu<<^, 



which holds for the given weights. □ 

This then gives us an expression for the matrix of the self-adjoint operator X*T 
with respect to our orthonormal basis of polynomials: 
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(4.17) = ^^iVf-". 

We next define a matrix ^ that will give coefficients of the functions Vk ° t 
expanded in terms of the monomials {wjljeNo- 



(4.18) VkOT = ^AjkVj. 

3=0 

The entries of A exist since we assumed that each Wfc o r is an element of V and 
therefore has an L^-convergent expansion in the monomials. We may or may not, 
however, be able to compute entries of the matrix A directly from this definition. 

Example 4.6. A nonafhne map: t{x) ~ + b 

This is perhaps the simplest example on a nonaffinc transformation. The matrix 
A which encodes t{x) = + 6 can be computed from the powers of t to satisfy 
Equation (|iT5)) : 



1 


b 


62 


63 

















1 


26 


362 
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We wish to think of A as the matrix representation for an operator on the 
weighted space Vw- With an abuse of notation, we will also refer to this operator 
as A. In general, it is not even certain that the operator A has dense domain. 
Throughout the remainder of this section, however, we will restrict our attention 
to the cases in which the domain of A contains the monomials, and hence the 
orthogonal polynomials. 

Proposition 4.7. Given A as in Equation {4^.1^ and such that A is the matrix 
representation of an operator (also denoted A) on L'^i/J.) with respect to the ONB 
{Pk}k<£Noj if the operator composition RA is densely defined, then 

RA^T 

on the domain of RA, i.e. RA is a restriction of T . 
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Proof. Suppose / e dom(i?). Then by Parseval we have 



oo 

Rf = ^'<Pj\I)LHp.)RPj 

oo 



with convergence in the sense. 

We have assumed that, as an operator, A has the monomial functions {vk}k 
and the polynomials {pfcjfceNo in its domain. We have also assumed that the product 
RA is well-defined on a dense subset of L^iiJ,)- We can then compute 



Tpk = VkOT 

oo 



OO 

= RApk. 

The last line above holds for every Apk that is in the domain of i?. □ 

Theorem 4.8. Given {X,fj,) a Borel measure space with X C M and given r 
a measurable map from X to itself, let N'^'^^ he the adjusted moment matrix for /i 
and N^^°'^ ' the corresponding moment matrix for ^ o . Then, if the matrix A 
from Equation \4-lS^ satisfies the hypotheses in Proposition \4.7\ 

(4.19) A*N^^''>A^ N^^'°^'">. 

Proof. We have RA = T. Find weights w such that both and R^, are 
densely defined on L^(/^). Let A^™'"' be the weighted moment matrix with entries 
^wjw i^il'^j) L^ii^) Using Equations (I4.14p and (|4.17|) . we find that 

w w 

= (RAYiRA) 
= A*RIRA 
^ A*Ni^^A. 

□ 



4.2. Approximation of A with finite matrices 

In this section, we explore whether we can perform computations with finite 
matrices which yield a finite approximation to the infinite matrix A defined in 
Equation (|4.18p . and thereby achieve an approximation of the moments for /io r^^ 
from the moments for fi. 

Let /i be a Borel measure on a set X C M such that the moments m/^"* = 
j^x^^^ dii{x) are finite for all orders. Let t : X ^ X he & measurable endomor- 
phism such that the moments with respect to /i o r^^ are also finite for all orders, 
and the powers of r are in the closed span V of the monomials in L'^ifJ.)- Let T 
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be the infinite matrix introduced in Section [4.11 with entries Tij ~ {pi\vj 
where {pijigN^ are the orthonormal polynomials in L'^ip) given by performing the 
Gram-Schmidt method on the monomials {w^l^gNo- Let R be the transformation 
taking pi to Vi, so the matrix entries are Ri^j = {pi\vj) i^-i^^^y 

Fix n. Let i?„ and Tn be (as in Section [3?2)) the (n + 1) x [n + 1) truncations 
of the R and T matrices, respectively. Let Hn be the closed linear span of the 
monomials {vq, vi, . . . , Vn} (which is also the closed linear span of {pj;}"^o) and let 
P„ be the orthogonal projection onto 7Y„ from V. Note that in the Dirac notation 
mentioned in Chapter [1] 

n 

(4.20) P^^Y,\P^){P^\■ 

i=0 

We have proved in Proposition 14.71 that the matrix A = R~^T, in the cases 
where this product of infinite matrices is well defined, gives the coefficients of Vj o r 
in terms of the monomials {wijigNo- We now wish to show that the finite matrix 
product R~^Tn provides an approximation of A, in the sense that it yields coef- 
ficients for the projection of the powers of t onto Tin, expanded in terms of the 
monomials. 

Lemma 4.9. For fixed n, 

n 

(4.21) V(i?,7ir„),,fc«j = P„(«fc o r). 



3=0 



Consequently, 



(4.22) lim y (i?~'T„),- fcw, = Vk o r, 

where convergence is in i^(/i). 

Proof. Given our fixed n, let k < n. In the computation below, recall that 
for j, k < n, wc have {Rn^)j,k = RJk ^^'^ {Tn)j,k = Tj^k- Also, since the orthogonal 
polynomials po, . . . ,pj are in the span of the monomials vq, ■ ■ - Vj, we know that the 
truncated matrix R^^ maps Vj to pj for each j < n. 



Y^iR-'T^hkV, = y5](i?„)-l(T„),,,«, 

n n 

= y o T)i2(^) y (i?,7^)j,fUj where both e,j< 

1=0 3=0 



^{pe\vkOT)L2f^^-j 



e=o 

= PnVk °T by Equation (|4.20p . 

□ 

Remark 4.10. It is important to realize here that even for j, k < n, 
(4.23) {Ru'Tnhk + {R-^T\k. 



40 



4. MOMENT MATRIX TRANSFORMATION: MEASURABLE MAPS 



We do, however, know that because the product R is well defined, if we fix j, k 

and let n ^ oo, we do have {R:^^Tn)j,k {R^^T)j^k- 

If we were able to conclude from the truncation result in Lemma 14.91 that 
lim„^oo X]j=o(^~^-^)i.*:^j = Vfc o T, where wc take the limit in we would 

have an alternate proof of Proposition 14.71 It other words, we would be able to 
write 

oo 

(4.24) Y.^R-^^)j,kVj = Vk o r. 

j=o 

This convergence, however, would require hypotheses about how the sequences 
(i?^^T„)j.fc converge to {R^^T)j^k as n — > oo. To illustrate this point, let us fix k 

and let a^"'' = {R^^T„)j^k and aj = {R^^T)j^k, so for each j, a^"^ — > aj as n oo. 



j=o j=0 j=0 j=0 

n n 



II ^ ^ lli^tp) II ^ 

3=0 3=0 



Using Lemma 231 the second term of the last line above can certainly be made 
arbitrarily small for large enough n. The first term, however, may not have that 
property. 

For each fc e No, the truncated matrix products R~^Tn produce asymptotic 
expansions for Vk o r by giving expansions of Pn{vk ° t) in terms of monomials Vj, 
even if there is no a priori known expansion of t^fe or in the monomials. We assume 
that such an expansion exists, that is, 

oo 

Vk o T — aj Vj . 
i=n 

We may, however, have no way of computing the actual coefficients. Lemma 14.91 
gives approximations to these coefficients which get better as n increases. The 
next question is how these finite approximations to the matrix A interact with the 
moment matrix transformation given in Equation (j4.19p from Tlieorem l4.8l In the 
case described in Subsection 13.11 where t is an affine map on M or C, the matrix 
A is upper triangular. With this added structure, it is readily demonstrated that 
the truncated matrix product yields exactly the truncation of the infinite matrix 
product, i.e. 

We also note that in this special case, we have Equation p.5p from |EST06j giving 
a concrete expression of the entries of A. We then can compute the exact moments 
with respect to ^ o using finite matrix computations. 

In the more general case, however, we may not have a construction giving us 
the entries in A, so we may need to use the finite matrix product An = R^^Tn- In 
this case, we find the triple product gives entries which are the inner products of 
projections of powers of the measurable map r. 

Given the map r, the entries in the adjusted moment matrix for /i o are 
given by the inner product 



4.2. APPROXIMATION OF A WITH FINITE MATRICES 41 

Corollary 4.11. Let An = i?,7^T„. Then 
(4.25) (I;7V(^)I„),,, = (P„(t;, o t)|P„(«, o T))i2(^) = (t;, o t|P„(«, o T))i2(^). 
Proof. This is a consequence of Lemma [4.91 

n n 

?i n 

= ^^{An)kA'"k\vi)L^{^){An)lJ 
k=0 1=0 

(n n \ 

k=0 1=0 I L'^(^) 

= {Pn{ViOT)\Pn{Vj Ot))l2(^^,) 
^ (Ui O t|P„(Uj O t))l2(^). 

In the last line, we use the property of projections that Pn — Pn — Pn- ^ 



CHAPTER 5 



The Kato-Friedrichs operator 

The moment matrix M^''^ is an infinite matrix, and while M'^^^^ may not be a 
weU defined operator on we may be able to view as an operator in some 

other sequence space. We use the techniques of Kato and Fricdrichs to turn A/'^' 
into a self-adjoint densely defined operator on a weighted l'^ space. Our main tool 
will be a quadratic form Qa/ which is defined from the moments of /i. We obtain the 
weighted l"^ space by finding a space in which the quadratic form Qm is closable. 
In the course of showing that Qm is closable, we introduce two key operators: F 
and its adjoint F* , which we will continue to study in Chapters [6] and [71 

5.1. The quadratic form Qm 

Let T) be the set of all finitely supported sequences indexed by Nq. The familiar 
expression in Equation (|1.12p for the infinite matrix- vector product Mc makes sense 
for every c ^T), but Mc may not be well defined for every sequence c. Even when 
it is well defined, Mc may not be in the same Hilbert space as c. We will address 
this technicality by changing the domain of M . 

Let M be the moment matrix of a positive Borel measure with finite moments, 
where M satisfies the positive semidefinite condition in Definition 12.11 (real case) 
or the PDC condition in Definition 12.31 (complex case). The quadratic form Qm is 
defined on infinite sequences with only finitely many nonzero components. Specifi- 
cally, given c € I? we have 

(5.1) Qm{c) = ^ ^ c7MijCj. 

i j 

Even though every operator on a Hilbert space determines a quadratic form, 
the converse is not necessarily true. Here, we look for conditions on the quadratic 
form Qm which guarantee that the quadratic form gives rise to an operator on a 
Hilbert space Ti. The work of Fricdrichs and Kato provides the correct conditions. 
Friedrichs (see jKatSOl Section 6.2.3]) proved that semibounded symmetric oper- 
ators in a Hilbert space Ti. have self-adjoint extensions in Ti. Later, Kato [Kat80[ 
Section 6.2.1-Theorem 2.1; Section 6.2.6-Theorem 2.23] extended Friedrichs's the- 
orem to quadratic forms. 

Theorem 5.1 (Kato). Let Q he a densely defined, closed, positive quadratic 
form. Then there exists a unique self- adjoint operator H on a Hilbert space Ti. such 
that the Hilbert space completion Tig of Q is equal to the completion of H , and for 
all c in the domain of Q, 

Q{c)^\\H^'\rn = \\c\\l,^. 
In particular, the domain of Q is equal to the domain of H^^^. 
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Our quadratic form Qm arising from a moment matrix may not be closed, but 
we can show that Qm is closable and then apply Kato's theorem to the closure. 
We change our notation briefly here in order to work with sequences of elements 
in D. A sequence in D will be denoted {c„}„gNoj ^-nd the i"^ component of {c„} is 
Cn(i). (We retain the notation c ~ {c{i)}i^-Mg for elements of V and its subsequent 
completion throughout the next two sections.) 

Definition 5.2 ( [KatSOj ). Let H be a Hilbert space in which V is dense. A 
quadratic form Q defined on 2? is closable if and only if 

c„ in TC 

and 

Q{Cn — Cm) ^ as TO, n ^ oo 

imply 

Q(cn) ^ as n — > oo. 

We next define the operator F to map finite sequences c e I? to polynomials 
(a.k.a generating functions) in L^(^). Given c e we define 

(5.2) Fcix) = f,{x) = J2 

ieNo 

Example 5.3. A quadratic form which is not closable in 

In order to construct a quadratic form which is not closable in we need a 
sequence {c„} C i'^ and a measure such that 

(1) Ei6No |cn(i)P ^ as n ^ GO 

(2) /k l/c™(a;) - /c„(a;)pd/i(a;) as m,n^ oo 

(3) /jj |/c„(a::)p d^(a;) 7^ as n -> oo. 

Let ^ = (5f, for 6 > 1, and let c„(i) = 5{n,i)h^^ (the Kronecker delta). Since 6 > 1, 
||c„||p = 6^^" ^ as 71 oo, so (1) is satisfied. Since 

/c„ {x) ~ 6 "a:^" for each n, 
and the integral in (2) simply evaluates /c„ and fc^ at x = b, we have 

/ \fc„M - fcAx)?dfi{x) = |1 - 1|2 = for aU m,n. 
Jr 

However, the integral in (3) is 1 for all n. Therefore the quadratic form associated 
with 6b, b > 1, is not closable in 

5.2. The closability of Qm 

In order to show that the quadratic form Qi\i is closable for a particular infinite 
Hankel matrix Af, we show that the operator F defined in Equation (|5.2p is closable 
with respect to a weighted space £'^{w). Recall Definition II .21 regarding closed and 
closable operators. The weighted P' spaces are defined by a sequence of weights 

e{w) ■.= [{c{i)}\Y^w,\c(S)f 

The norm on i?^(w) is 

l|c||^2(^) = ^ w,|c(^)|^ 

i6No 
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and the inner product is given by 

{c\d)i2(^) = ^ Wic{i)d{i). 

We will generally take the weights to all be strictly greater than zero. 

In the next lemma, we use an argument mirroring the proof of Lemma 14.41 to 
produce a set of weights w which guarantee that F will be closable in the weighted 
space (^{w), and then in Theorem 15.51 wc show that the quadratic form Qm is 
closable in P{w). 

Consider _F as a map with domain in P('w) for some weights w = {wi} nzfi^. 
We then see thcit the adjoint F"^ of F is given by 

(5.3) (F;/)fc = — /xV(a:)dMx). 

Wk Jr 

To verify this, wc must show that F^ satisfies 

(5.4) {Fc\f)m^.) = (c|F:/),2(„) 

for every / such that F^f £ i'^{w). 
For each c £ T), we have 



(5.5) 



{Fc\f)L^,.^= [ 5] c(fc)a:V(x)dMx) 

I xV(x)d/i(x). 

',.^^T J«. 



We now multiply and divide by Wk'. 

Jr 



(5.6) 



keNo 



where we have shown the right had side of Equation (|5.3|) satisfies the definition of 
the adjoint F^. 

Lemma 5.4. Suppose |a;|'^' d/.t(x) < oo for each fc G No. There exist weights 
{wi}i£N„ such that F : Piw) i^(/^) ( t5.^|) is closable. 

Proof. The operator F is closable if and only if dom(F*) is dense in i^(/i). 
(See, for example, [Con90[ Proposition 1.6, Chapter 10].) 
Now, F^f € P{w) precisely when 

(5.7) 5^«;fc|(F:/)fcp<oo, 

fceNo 

and (|5.7p is true if and only if 

2 



(5.8) E— / ^''fi^)Mx) 

u ^"k Jr 



< oo. 



k 

This gives us one sufficient criterion for F to be closable. We see from Equation 
(j5.8|) that the monomial functions {wjjjgNo ^re in the domain of F^ if and only if 

(5.9) T.V-\I ^■'■+'d/i(x)' = ^— |M,,fc|2<oo 

^ Wk\ Jr ^ Wk 
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for all j e No. 

Rather than looking at this criterion (j5.9p . however, we will find weights that 
put the essentially bounded functions in the domain of F^. By hypothesis, /i 
is a finite measure so L°°{^) C L^{^). Set Aik ■= Is;]*^ d/i(a;), and suppose 
/ e i°°(At)- Then 



(5.10) / x^fix)dti{x) < \\f\\L^^^)Mk. 

Jr 

Therefore, if {Ml/wk} € then 

(5.11) L°°(m) C dom(^^;), 
and we explicitly define a choice of weights Wk ■ 

(5.12) Wk := {I + k^)Ml 

Since L°°{fi) is dense in L^{n), we now know that dom(i^^) is dense in £'^(w) where 
w — {wk}ketio is defined in (j5.12p . Therefore F : £'^{w) ^^(m) is closable. □ 

Theorem 5.5. Given fi a Borel measure on R with finite moments of all orders 
and M — Af'^\ there are weights {ti;(i)}igNo such that {Qm,'D,£'^{w)) is closable. 

Proof. Choose weights w = {wijigNo a-s in Equation (|5.12p so that F is 
closable. Suppose c„ — > and QM{cn " Cm) as m, n — > cx) as in Definition 

15.21 Our goal is to show that QM{cn) ^ 0. 

Denote the polynomial Fc„ by /„ for each n € Nq. Because Qm on T) and Qm 
defined on the polynomials satisfy 

Qm{F{c)) = Qm(c), 

the condition QMic-n — Cm) ~* implies that {fn}neKa is a Cauchy sequence in 
L^ifJ.)- Therefore, there exists g E i^(/i) such that /„ — > g. 

Wc will now use the condition that c„ in £'^{w). By Lemma 15.41 F : 
P{w) — > i^(A*) is closable, so its closure F has a closed graph. Since c„ ^ 0, 
F{cn) = F{cn) g, and the graph of F is closed, g must be S L^il^)- Q 

As a result of Theorem l5.51 we can apply Theorem l5.1l (Kato's Theorem) to the 
closure of the quadratic form Qm- We denote by Tiq the Hilbert space completion 
with respect to the quadratic form Qm- By Kato's theorem, there exists a self- 
adjoint operator H on (."^{w) such that the Hilbert space completion is also Tig. 
This gives 

(5.13) Ml,^ = Qm{c) = \\H''^c\\l^^~, 

for all c e Ti-Q- Moreover, the domain of H^^^ is equal to the domain of Qm- It 
should be noted that in general, the domain of is a subset of the domain of H^^^ 
and may or may not contain T). 

5.3. A factorization of the Kato-Friedrichs operator 

Given F from £^(w) to L'^ifJ,), suppose the weights w have been selected such 
that F is closable, i.e. F^ is densely defined in i^(/x). 
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Let Qm be the quadratic form given by the i^(//)-norni: 

(5.14) QM{f)= I \f?d^i. 
Then for every c G 2?, we have 

QaAFc) = ! \Fc\''d^Ji 

JR i J 

(5.15) = ^^ciMi^jCj = Qm{c). 

i 3 

Proposition 5.6. F'^^F is a self-adjoint densely defined operator on £'^{w), and 
F^F is the Kato operator H . 

Proof. The domain of F contains V, and the domain of F*^ contains the 
polynomials. Since F maps finite sequences c € I? to polynomials Fc, we sec that 
V is contained in the domain of F^F. In order to prove that F^F is the Kato 
operator corresponding to the quadratic form Qm, we must show that it satisfies 

Qm{c)^\\{F:F)^/'c\\,.^^^ 

for all ceV. 

First, note that for any c E V, 

{c\F:Fc)e^^^ ^ {{F:Ff'^c\{F:,Ff'^c)^,^^^ 
= \\iF:Ff/'c\\%^^y 

We also note that 

(c|f;Fc),2(^) = {Fc\Fc)e2^^) 
= \\Fc\\l 

= Qm{Fc) = Qm{c). 

Therefore, for c G V, 

Qm{c) - {c\F*^Fc),.(^) = ll(^:^^)''''c||f.(„). 

By uniqueness of the Kato operator for the closable quadratic form Qm, we have 
F^F = H. □ 

This proposition makes the usually mysterious Kato operator much more con- 
crete. In particular, as a result of this proposition, the domain of H includes V. 

5.4. Kato connection to A matrix 

Given a moment matrix M^''' and a measurable map t, we can use Kato theory 
as another way to describe a matrix A such that, under appropriate convergence 
criteria. 

This is a different approach to that used in Section [4?Tl We prove that under certain 
hypotheses on the measure fj, and the map r, we can find an operator A which is 
an isometry between the Kato completion Hilbert spaces for the quadratic forms 
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Qm and Q'j^j, where M = Af(^) and M' = A/(^°^ We demonstrate that this 
isonictry A wiU intertwine the Kato operators for Qi^i and Qa/s a-nd therefore the 
corresponding matrix representation A will intertwine the moment matrices M and 
M' as required. 

Let (X,^) be a measure space with X C M. For each k e No, let Vk be the 
monomial function on X: Vk{x) = . The polynomials may not be dense in L'^{n), 
so denote by V the closed span of the polynomials in We take /i to be a 

measure such that {wfe}fegNo ^ -^^(m) ^ ^- Let M ~ M^^^ be the moment matrix 
for n. By Section l5.ll there exist weights w = {wijigNo such that the map F is 
closable on the weighted space Thus the quadratic form Qm is closable, and 

we can define the Kato operator H on £'^{w). Note that the domains of Qm, F, and 
F^ all contain the set of finite sequences V. 

Let r be a measurable endomorphism on X such that the functions Vk°T — 
are all in L^{fJ,) n V. Let M' ~ A/'^""^ ' be the moment matrix for the measure 
fi o T~-^. Define the map F' : V L^U^- ° t^^) to be the analog of F, i.e. the map 
which takes c & V to the polynomial fc = J2i '^i'^^i ^ o t~^). Using the same 
reasoning as that used in Section 15.11 there exist weights such that F' is closable, 
hence the quadratic form Qm' is closable, and we have H' , the Kato operator on 
the corresponding weighted space. 

Let w = {wfejfcGNo be weights that allow for both F and F' to be closable on 
£'^{w), hence the Kato operators H and H' are both densely defined self-adjoint 
operators on £'^{w). Denote the Kato completion Hilbert spaces for the quadratic 
forms Qm and Qm' by Hq and Hq' , respectively. 

Notice that S i^(/io t~^) if and only if Jj^ |0o r^d/x < oo, which holds if and 
only if (/) o r £ L^{fi). In fact, there is a natural isometry a between o t~^) 
and ^^(/x) given by 

a(0) — (j) o T. 

Remark 5.7. The map F which we have defined above actually takes on several 
realizations in this section. Simply put, it is the map that takes finite sequences 
to the polynomial with these coefficients. But we can think of F as operating 
on the Kato spaces Tiq, Hq' or on the weighted space We can also take 

the codomain of F to be L^{^i) or L^{^ o t~^). In order to clarify these different 
realizations, we will denote the operator F as follows: 









Fq' 


T^Q' - 




F 






F' 


£\w) 





Moreover, we also observe that the operator T (defined below) which maps standard 
basis elements to can be expressed in terms of the F' and the isometry a. 
Given c G P, 

Tc = (5]c,r^) 
= aF'c. 
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We will show below that the operators Fq and Fqi are isometries, and we have 
already proved that for appropriately chosen weights, the operators F and F' are 
closable on £'^{w). 

Lemma 5.8. The Kato space Hq is isometric to V Q Lp'{^), andT-Lqi is iso- 
metric to a subspace of L^(fi o t~^). 

Proof. The norm on the Hilbcrt space Hq is given by 

Mhq = Qa/(c) 

when c £ V, and by the definition of the Kato operator, we also have 

\\c\\n^ = \\H'/'c\U2^^y 
Let Fq be the map taking c E V to fc = '^i'^i in 7^ C L^{fi). Then 
II^^Qcll^ = WfcWl 



2 



i 

Qm{c) - ||C||2 . 



Thus Fq is an isometry from T-Lq to V . An identical argument proves that 
the corresponding Fqi map from Tiqi to o r~^) is also an isometry onto the 
subspace V' of o r^^) spanned by the polynomials. □ 

Proposition 5.9. The Kato Hilbert spaces T-Lq and Ti.Qi are isometric. 

Proof. We use the notation = if there exists an isometry from the Hilbert 
space Ti. onto the Hilbert space Q. Using Lemma 15.81 and the map a, we have the 
following isometries: 

tcq- ^v = nQ. 

□ 

If we denote the isometry in Proposition 15.91 bv A, then 

A^FQ^aFQ,. 

The isometry property \\c\\q' = j| j4c||q implies the quadratic forms satisfy Qm(Ac) = 
Qm'[c) for all cEV. 

Given {eijigN^ the standard orthonormal basis for we define the standard 
orthonormal basis {ef }igNp = for t'^{w). Given F from (."^{w) to 

suppose the weights w have been selected such that F is closable, i.e. F^ is densely 
defined in L^{^)- Recall from Lemma [5.61 that F^F is equal to the Kato operator 

The operator A has domain containing T> C TYg, but we will need to consider 
the analog of A as an operator on £'^{w). In order to do this, we define the following 



50 



5. THE KATO-FRIEDRICHS OPERATOR 



maps: 



Hq 



Each of jw,j'w ^cts like the identity map on T>, mapping in ^^(w) to in the 
Kato spaces Hq, TLqi, respectively. Then, we see 



and 



F' = Fq,j'^ 



We then wish to define the operator on £'^{w): 

Remark 5.10. Defining A in this fashion requires the composition of opera- 
tors to be defined on T) and also requires jw be invertible. We will assume these 
properties through the remainder of this section, but we note that they may not 
hold in general. Thus, this approach requires convergence hypotheses to define A, 
just as were needed in Proposition |4? 



Given the maps a and F' defined above, we define T to be the composition 
T = aF' : fiw) -> L^{fi), 
where T maps basis element to the function ^^t' = -^=Vj o r. 

Wc now wish to write T with respect to the isometry A and the map F ~ FqJ^ 
We will need to make use of the operator A, as shown in the following diagram. 



Hq' 



Hq 



Lemma 5.11. Given the operator T defined above, 



FA. 



Proof. Recall from the definition of A that aFg = FqA. We expand T using 
the diagram above. 

= c^f^jL 

= FqAJ'^ 

= Fj-^AjI^FA. 



□ 



Lemma 5.12. There exist weights such that T is a closable operator on the 
weighted space £'^{w). 
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Proof. We have proven that there exist weights w = {wijigNo aheady that 
make the operators F and F' closable, and in particular, these weights ensure that 
the bounded functions are in the domains of both F^ and (-F!^,)*. Since a is an 
isometry, these weights also ensure that has dense domain. □ 

Lemma 5.13. Let w be weights such that F and T are both closable operators 
on i^{w). If the operator T^T is self-adjoint and densely defined on £'^{w), then 
T*T^H'. 



Proof. This proof repeats the argument from Lemma l5.6l 



= —^={vi\Vj)L^^,or-^) 

Therefore, we have 

Qm'(c) = {c\M'c)m^) = ll(T*T)i/2c||2,^^^ 

for all c G V, which proves T*^T is exactly the Kato operator H' for Qm'. This 
proves T*,T is self- adjoint and densely defined on £'^{w). □ 

Let A be the operator on i'^{w), so T = FA. 

Theorem 5.14. The operator A intertwines the Kato operators H and H' , and 
thus in matrix form intertwines the moment matrices Af^^^ and M^'^°'^ •* . 

Proof. If the product T = FA is defined on P, we have 

T^T = {FA)* FA = A*{F;^F)A, 

which by Lemmas 15.61 and 15.131 implies that 

H' = A* HA. 

If we form the matrices for A, H, and H' with respect to the orthonormal basis 
{e™}igNo, we have 

A*M^^''>A = M^^'°^''l 

□ 

Using the Kato approach we can therefore show that, under the appropriate 
conditions on t and the measures /i and /i o r, we can find a matrix A such that 
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5.5. Examples 

The following examples illustrate the situations in which infinite matrices do 
not have direct realizations as i'^ operators, but can be realized as operators on a 
certain weighted Hilbert space. The simplest such case arises when the measure /i 
is a Dirac mass at the point 6=1, which we observed earlier in Example 11.81 

Example 5.15. The Kato operator for ^ = 5i. 

Let /i = (5i be the Dirac measure on R with point mass at 1. As we stated in 
Example [LH the moment matrix M = M^^^ will be Mij = 1 for all i, j € Nq. The 
associated quadratic form Qm is given on D by 

I ^ |2 

Qm{c) = X! '^'^ • 

fcSNo 

If we attempt to define the operator F* on the constant function f{x) without 
weights, we have for each i, 

{F*f\= [ f{x)x^dS,{x)^f{l). 

The weights are necessary in order to make any function with /(I) ^ 1 be in the 
domain of F* . 

Select weights w = {wi} C such that 



E- 

IDA 



< oo. 

Wi 

Then F^ and Qm are closable, and the resulting Kato-Friedrichs operator then 
satisfies Equation (|5.13p . Given c is in the domain of H^, we have 

j k 3 

for all j e Nq. Thus, is defined by 



ID: ^ ^ 



Cfc. 

In other words, is a rank-one operator with range equal to the span of the 
vector 

> LlilQ Wi W2 1 



Example 5.16. The Laplace transform and the operator F* associated with 
the measure da; on IR+ 

We first observed in Examples 11.91 and 12.121 that the moment matrix entries 

M^,J = {i+jy. 

grow quickly as i and j increase. Therefore, we again need weights to make F a 
closable operator. 

Let w = {wijigNo be weights given by 
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We claim that with these weights, F is closable. 




[(2^)!]2 



for each fixed j G Nq. Using the ratio test, we see that the series converges for each 
j so by Equation (j5.9p we have the monomials in the domain of F*,, which makes 
F closable on £'^{w). 

As a side note, we observe that F^ has a connection, under appropriate condi- 
tions on the function /, to the Laplace transform C: 




where 




e 



— sx 



for s e M+. 







CHAPTER 6 



The integral operator of a moment matrix 

In this chapter, we will first show how the Hilbert matrix M, which is the 
moment matrix for Lebesgue measure on [0,1], can be associated to an integral 
operator. It turns out that both the operator associated to AI and the integral 
operator are bounded. In fact, we use properties of Hardy space functions and the 
polar decomposition to show M and its integral operator are unitarily equivalent. 
The bounded operator properties of the Hilbert matrix are well-known. See [Hal67l 
IWid66| for more details. 

We will generalize Widom's technique, using a weighted Hilbert space when 
necessary, to find an integral operator on L^{n) associated with the moment matrix 
M^^'\ where /i is a positive Borel measure supported on [—1,1]. In the more general 
setting, the integral operator may not be bounded. In Chapter [7] we will examine 
conditions under which the integral operator is bounded. 

When we refer to the Hilbert matrix, where yU is Lebesgue measure restricted 
to [0, 1), we will use the symbol M; for any other measure ^ supported in [—1, 1], 
the moment matrix will be denoted 



6.1. The Hilbert matrix 

The Hilbert matrix M is the moment matrix for Lebesgue measure restricted 
to [0, 1). M has («, j)*^ entry given by 

1 

(6.1) Mi , := = / x'+^ dx. 

It is well-known |Hal67l IWid66| that the Hilbert matrix defines a bounded oper- 
ator on and the operator norm of M as an operator from £^ to is tt. 

Using our definition of 2? to be the set of sequences with only finitely many 
nonzero coordinates, recall the operator F : 2? — > ■p[0, 1] given by 

Fc{x) = Cna;". 

neNo 

The operator F defines a correspondence between the finite sequences 2? and the 
polynomials on [0,1], and this correspondence extends by Stone- Weierstrass and 
the Riesz-Fisher theorem to an operator from £^ to L'^[0, 1]. Wc will often call the 
image of a sequence c a generating function Fc = fc. 

Suppose / S L^(0,1) and c € i'^. By a Fubini argument and the Cauchy- 
Schwarz inequality, we can switch the sum and the integral in (/|i^c)i2(o^i) to 
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define the adjoint operator F*: 

{f\Fc)L2 = 



/ /(a;)Vc„x"dx 

Vc„ / f{x)x"dx 
Jo 



n=0 

{F*f\c 



where we have the adjoint now defined by 



Jo 



Furthermore, if / G L^{0, 1), then {(i^*/)n}«eNo S 

Lemma 6.1. Let M be the operator for the Hilbert matrix, and let F and F* 
be the associated operators defined previously. Then F* F ~ M . 

Proof. Given c e the matrix product Mc is well defined since M is 
a bounded operator on We can therefore write {Mc)i ~ Sjlo -^(*' ~ 

E°° n ■,'^-•'11 , where convergence is in the l"^ norm. Then 



{F*Fc)., = [ x'{Fc){x)dx 
Jo 
„i 00 

= / x^'S^CjX^ dx 
Jo 



00 „i 



J=0 

x'''^-' dx 



00 

The exchange of the sum and integral above follows from Fubini since c G and 
i+j+i ) ^ '■^^ ^ sequence in j). □ 

We now turn to a relation between the operator M and an integral operator 
K. This relationship is shown in a more general setting in |Wid66| . and we will 
also discuss a generalized version in Section 16.21 

Theorem 6.2 f [Wid66| . Lemma 3.1). Suppose M = F*F is the operator 
representing the Hilbert matrix. The self-adjoint operator K = FF* on Li^[Q, 1] is 
an integral operator with kernel 

1 

k{x,y) 



l-xy 

Moreover, K and M satisfy the relation on P' : 
(6.2) KF = FM. 
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Proof. Let / e ^^[0, 1] and fix x e (0, 1). Then 

OO 

{FF*f){x) = J2(P*f)jX^ 



,. -n "'0 



3=0 
r* 1 ^ 

1-xy' 

Wc again use Fubini to exchange the sum and integral above, so FF* is the desired 
integral operator. Combining this equation with Lemma l6.ll we find 

KF ^ FF*F = FM. 

□ 

We next develop tools which we will use in the next section to study the spec- 
trum of the integral operator K. We will demonstrate the relationship between the 
integral transform and the Hardy space IHI2 on the open unit disc D in C. 

The Hardy space 1HI2 consists of functions which are analytic on D. The Hardy 
space norm of a function F{z) with expansion F{z) = X]^o ^n^" is given by 



The inner product on IHI2 is exactly the P inner product of the expansion coefficients 
of two Hardy functions. 

Given z € D. we let kz{x) ~ (1 — zx)^^ . Then K denotes the integral operator 
with kernel , which is the extension of K defined in Theorem 16.21 to the complex 
unit disc D. Given g e L^[0, 1], we define 

{Kg){z) = [ f^^l dx = (fcz|g)L2[o.i]- 
Jo ^ — zx 

Note that as a function of z G D, is analytic for x G [—1, 1]: 

nGNo 

Next, we prove that the range of K is contained in EI2. In particular, we show 
that if / = Kg^ with g G i^[0, 1], then / can be extended to an analytic function 
F{z) on the open unit disc D. 



Lemma 6.3. Suppose f G £^[0, 1]. Then 



dx G 
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Proof. Given x e [0, 1] and \z\ < 1, we have 

„1 oo 

F{z)= / /(x) Va;"z"da;. 

•^0 n=0 



Next, observe that (/^ \f{x)\x^ dx) is in i'^ as a sequence in n, since by Cauchy- 
Schwarz 



/I|l2[0,1] 

2n+l 

This aUows us to use Fubini's theorem to write 



\fix)\x-dx = (|/||a;")i.[o,i] < |i/|lL^[o,i] lla;"|U.[o,i] = 
3ini's theorem to write 

n=0 ^-^0 / 

and also shows that F e EII2. □ 

The adjoint operator K* to our extended integral operator JsT is a map from 
H2 to i^[0, 1]. Let (f> GM2 have expansion 0(z) = J2'^=o Cn^:". Then, 



= / f{x)'^CnX'^dx 



n=0 

= (/I0>L2[O4]- 

Wc use Fubini's Theorem above to switch the sum and the integral, and note 
that in the last line, (j) is restricted to the real interval [0,1]. Thus, wc find that 
K*(j) is the restriction of cj) to [0, 1]. We now provide an argument using the polar 
decomposition and Hardy spaces to show that the operator F is in fact a unitary 
operator, so the operators M and K have the same spectrum. Using the definitions 
of F and F* , we note that the polynomials are in the domain of F* . This proves 
that F* has a dense domain, and thus F is closablc. Wc therefore have by a result 
of von Neumann that both F*F and FF* are self-adjoint positive operators with 
dense domains. 

By polar decomposition, given the closable operator F : £^ — ^ L^[0, 1], there is 
a partial isometry U : ^ 1] such that 

F = U{F*Ffl'^ = {FF*y^^U. 

Because C/ is a partial isometry, it is an isometry on the orthogonal complement 
of the kernel of F. Our goal now is to show that U is, in fact, a unitary operator, 
i.e. the kernels of F and F* are trivial. This will give us a unitary equivalence 
between the operators M and K. 

Given c £ let Fc = 0, that is, J2'^=o Cna^" = for a.e. x in [0, 1]. But, by 
the definition of Hardy functions, X^^o Cn^:" is the restriction to [0, 1] of an IHI2 
function 

00 

n=0 

Since F is analytic and is zero on a set with accumulation points, we have F{z) = 
on D. Therefore Cn — for all n, which gives ker(F) = {0}. 
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For a function / G L^[0, 1], assume that F*f = 0. Thus, /q /(a;).T" dx = for 
ah n G No- This means that / is orthogonal to the monomials, and hence to the 
polynomials. Since the polynomials are dense in L^[0, 1] by Stone- Weierstrass, we 
have that f{x) = a.e.x. Therefore, ker(i^) = {0}, and we can now conclude that 
the operator [/ is a unitary operator from £^ to L^[0, 1]. 

Given our definition of U from the polar decomposition, we have U{F* F)^^"^ = 
{FF*y/^U. This gives 

U{F*FY'^U* ^ {FF*y/^, 
and squaring both sides gives 

U{F*F)U* = UMU* = FF* = K. 

6.2. Integral operator for a measure supported on [1, 1] 

In this section we study measures /i with finite moments of all orders on M. We 
generalize the results found in Section [6. II for the Hilbert matrix. Given a measure 
/X, we see that the moment matrix Af may be realized in two ways: 

(1) as an operator H having dense domain in l"^ or a weighted space (."^{w) and 
acting via formal matrix multiplication by M^^^; and 

(2) as an integral operator K in the L^(/x)-space of all square integrable functions. 
While the resulting duality is also true in W^, for any value of d, for clarity we will 
present the details here just for d = 1. The reader will be able to generalize to 
d > 1. We will also assume our measures to be supported in [—1, 1]. 

Given the moment matrix Af*-^-*, define a quadratic form Qm as in Definition 
15.11 on the dense space T) consisting of finite sequences. Recall from Section 15.11 
that if the monomial functions {wfej^gNg are in H L^(/x), there exist weights 

w = {wijiGMo such that Qm is a closable quadratic form. (Note: in some cases 
such as the Hilbert matrix, the weights are not required.) Moreover, the operator 
F : £'^{w) L'^ilA which takes c G V to the generating function (a polynomial) 
/c(x) = J2j ^j^'^ is also made closable by the same weights. Then, as we showed 
previously, 

Qm(c)= / \fc{x)\^ dfi{x) 



for all c £ v. By Kato's theory, Q^/ has a corresponding closable linear operator 
H on £'^(w) with dense domain satisfying Qi\i{c) = \\H^^'^c\\'^2(^^^y In particular, the 
Hilbert space completion Tig with respect to the quadratic form Qm is the same 
as the completion with respect to the operator H. (See Remark 16.71 for another 
viewpoint on this space Hq.) We also showed in Lemma [5?^ that the Kato operator 
H is given by F^F, which shows that the domain of H contains T). Thus, for c G P, 
we have Qa/(c) = (c|i7c)£2(„,). 

We will show below that the operator FF* on i^(/i) is equal on its domain to 
the integral operator with kernel 

fc(.,,) = E^- 

Note that if weights are not required, the kernel above reduces to the same 
kernel which represented the Hilbert matrix in the previous section. The integral 
operators corresponding to moment matrices M'^'^^ for different measures ^ might 
have the same kernel, but they will act on different Hilbert spaces L'^{ij)- 
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We now can state the association between the Kato operator F*,F for a moment 
matrix M^f^\ the operator K = F^, and an integral operator on L'^{^jl). As usual, 
we use the notation Vj{x) = for the monomials in L'^di). 

Proposition 6.4. Let ji be a Borel measure on M with support contained in 
[— 1, 1] and Vj G L^{^J') for all j S Nq. Let F be the closure of the operator in i f <5. ^|) 
which sends a sequence c G T) to the polynomial fc G L'^(^). Let w — {wijigNo 
weights such that F^ has dense domain as an operator into £'^(w), i.e. via Equation 
(15. 9\) we assume that 



(6.3) V — |Af,+j2 < oo for all k e Nq. 

j<£No ^ 

Then the operator K = FF^ is a self-adjoint operator with dense domain. On its 
domain, K is an integral operator with kernel 

Proof. For the weights w, F is closable and the domain of F^ is given by 

(6.4) U€L\^,)\Y,w,\{F*^^),\' <^\. 



We note in particular that the monomials {uj}jeNo ^-re contained in the domain 
of F^ . The domain of FF*^ is the set of all functions in the domain of F^ such 
that F^{(f>) is in the domain of F. This set is dense by von Neumann's polar 
decomposition but may not contain all of the polynomials. 

Let (j) S dom(FF*,). Then a Lcbesgue dominated convergence argument gives 

iFF:)^ix) = J2 {F*<P)kVk{x) 
feeNo 

fceNo 



(^■5) JESl fct^„ Wk Jsuppip) 



I 0(y)^i:f^dMy) 



1PP(M) k&% 

The requirement that /j, is supported on [—1, 1] ensures that the sum above converges 
for all X, y. □ 

The polar decomposition of F is 

(6.6) F^UiF^F)'^^, 
where ?7 is a partial isometry. We therefore also have 

(6.7) F^{FF*Y'^U. 

The partial isometry U : l"^ L'^it') is the same in both (|6.6p and (|6.7p . 

Corollary 6.5. The two operators H = F^F and K = FF^ have the same 
spectrum, apart from the point 0. 
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Proof. A rearrangement of (|6.6p and (j6.7|) yields 
(6.8) UHU^K*. 

□ 

Example 6.6. Convex combinations of Dirac masses. 

Let fj, = ctiSxi, where ai > 0,^27=1 '^i ~ ^ ^^'^ ^i^^ ~ l,2,...,n are 

distinct real numbers. Recall from Example 1 5 . 1 51 that weights are required for any 
Dirac measure in order to make F closable. The same holds for convex combinations 
of these measures. It turns out, however, that the integral operator is not so difficult 
to describe. We show here that has dimension n, so the corresponding integral 

operator must have a representation as a finite matrix. 

Suppose the polynomial fc{x) = X^I-Lo ^kX^ is equal to the zero vector in L^{^). 
In other words, 



k=0 

Using our usual notation, fc is the image of a sequence c € under the map Fq 
from the Hilbcrt space Hq. Wc recall from Lemma [5.81 that Fq is an isomctry, 
hence 

\\fcfm^^=QM{c) = \\cf^^ 

where AI = is the moment matrix for /i. The matrix A/ is of the form 

n 



1=1 

Therefore, if /c = in L^{n), we have 

/ I ^ Ckx'' d/i(x) = ^ cJCkMj^k 



m n 
j,k=l i=l 

= XI "'I XI '^'^^i' 
i=l fc=0 

= 0. 

Since each > and X]r=i '^i ~ the above sum is zero if and only if 



k=0 



for alH = 1, 2, . . . , n. Therefore, the distinct real numbers roots of 

any polynomial which is equal to zero in L^(/i). 
Let p{x) be the degree-n polynomial 



P{x) = J|(x- a;^). 
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Using the Euclidean algorithm, any polynomial q can be written in the form q = 
ap + r, where a, r are polynomials and r has degree less than n. In the equivalence 
classes of L'^{ijl). wc therefore have that q = r. Hence, the dimension of L'^{ij.) is n. 



Remark 6.7 ( |j6ooj ). We can describe another way to generate the Hilbert 
space completion Tig of a quadratic form Q. So far we have looked at two ways 
to generate isomorphic Hilbert spaces TCq. The first is to complete the finite se- 
quences V with respect to the quadratic form Qm. The second is the closure of 
the polynomials in L'^{fi), where /i is the measure whose moments are the entries 
of M. Here, we briefly outline a third method to generate the same Hilbert space 
which consists of distributions. For now, wc will restrict ourselves to the concrete 
case M = the Hilbert matrix. 

As in our first method, we start with the positive semidefinite quadratic form 
generated by the moment matrix Af : 

Qm ■■ [0, 1) X [0, 1) ^ M. 

We use M to define linear functionals. For each x in [0, 1), there is a linear functional 
Vx defined by 

(6.9) Vx{-)^QMi-,x). 

(The analog to Vx is a sequence in V with zeros in every position except one, where 
there is a 1.) We use the VxS to build more linear functionals with finite sums (here, 
the analog is V): 

(6.10) «(•)= Qm{-,x). 

X finite X finite 

Finally, there is an inner product in the space of linear functionals of the form 

(6.11) Y c^'^^' Y CyWy^ = ^^c^CyQA/(a;,y). 

X finite y finite x y 

Now, we complete the inner product space to the Hilbert space Ji. The Hilbert 
space 7i is a reproducing kernel Hilbert space. It is fairly easy to check that for Vx 
as in (|6.9p and u as in (|6.10[) . (vx,u)j^ = u{x). 

Recall that the Hilbert matrix is related to the integral kernel (1 — xy)^^. We 
will next explain how this kernel appears in the Hilbert space of distributions TC. 

The tensor product of two distributions ui and U2 is straightforward; we apply 
the linear functional ui U2 to the pair of functions (/, g) by 

(6.12) {uT ® U2){f{x),g{y)) - W{f)u2{9). 

Then, we can extend this definition to functions h{x,y) because the span of func- 
tions of the form f{x)g{y) is dense in i^[0, 1]. 

One important example of a tensor product of distributions which is related 
to Widom's theorem is the following. Suppose 4>i and 02 are locally integrable 
functions on [0, 1), and ui = 4>idx and U2 = 02 dx. Then 
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The Hilbert space Ti, is spanned by distributions u on such that 

{u®u) \ — ^ — I < oo. 
\l~xy j 

It can be shown that the Dirac mass at and aU its derivatives belong to 7i; 
moreover, derivatives of the Dirac mass at (when scaled appropriately) form an 
orthonormal basis for Ti, with respect to the inner product (|6.1ip . The specific ONB 
is 

ir(") n 1 o 

Un = 1— 0^ ',71 = 0, 1,2, ... . 

n! 



CHAPTER 7 



Boundedness and spectral properties 



In cases such as the Hilbert matrix from Section 16. li the operators F and 
F* (which might be weighted or unweighted) are bounded operators, although in 
general they will be unbounded densely defined operators. In the first part of 
this chapter, we give sufficient conditions such that these operators are bounded 
operators, and hence the Kato operator F*F is also bounded. In Section [7751 we 
use the theory of projection- valued measures described in Section 17.21 to analyze 
connections among the spectrum of the Kato operator for a moment matrix M^f^\ 
the measure /i itself, and the associated integral operator. We will demonstrate 
these connections via the examples in Section [7. 51 We define the rank of a measure 
in Section 17.41 and demonstrate with examples which are convex combinations of 
Dirac measures. 



7.1. Bounded Kato operators 

Given a measure n with finite moments of all orders, denote its moment matrix 
by M = M'^^^ . We will use the generating function (x) having the even moments 
Mk^k as coefficients to express conditions which ensure the operator F^ on the 
weighted space P{w) is bounded. 

Proposition 7.1. Let fi be a measure with compact support onR with moments 
of all orders. Let the generating Junction 

oo 
fc=0 

have a finite and positive radius of convergence R. Select t such that t > R and 
J < R. If we choose weights Wk ■= t^ , then F^ is a bounded operator with 



I ^,!, 1 1 OB < Gu 



Proof. Note that since j < R, then j is within the radius of convergence, 
hence 

k=0 

We take Wk ■= t^ for our weights. 
Now, 

(7.1) (F;0)fc = — /xV(a;)d/^(x),fc € No. 

Wk J 
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So 



and 



(7.2) 



< M 



2k\ 



\F 



wV\\P(w) 



^ Wl 

fc=0 

oo _j 

< E -^^'^^ 



2 



fe=0 



GJ- 



Therefore, -F^ is a bounded operator with 



|-^lull£2(tu)_L2(^) 



sup ||i^*0||^2(„) 

11011 = 1 



<GJ- 



Using the properties of the adjoints of bounded operators, we also can conclude 
that F, F^F and FF^ are bounded operators with norms given by 



\f:\\1=\\f\\L^\\f*f 



wWop 



I op 



w-^ w\\op 



= \\F„F„ 



□ 



Lemma 7.2. Suppose fi is a measure on R with compact support and finite 
moments of all orders. Suppose in addition that supp{^) C [—1,1], with dfi{x) — 
B{x)dx, where B is a nonnegative bounded function. Then the operators F, F* , 
F* F and FF* are bounded, and 



\F*\\l, 



\F\ 



\F*F\\ 



\FF*\\op < 7t\\B\ 



lop 11^ ^ nop 

Proof. By }Wid66| (a straightforward generalization of Theorem 16. 2p . we 
know FF* is an integral operator of the form 



1 

-1 1 - 2;y 



dy. 



_i 1 - 

Let K be the positive integral operator defined in Section 16.11 which is equivalent 
to the moment matrix for Lebesgue measure. Then 

'■^ W)B{x)<P{y)B{y) 



• FF* 



l~xy 



■ dydx 



op\ 



\(l){x)B{x)\^dx. 



We recall that the operator norm of K is tt. In turn, the last expression is bounded 
above by 



ttWBW 



\cl>{x)\'B{x)dx = 7T\\B\\oo\mmf.y 
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Therefore, ||^^F*||op < ^||S||oo, and ||^^|| = < v/^M 



□ 



Theorem 7.3. Suppose there exists a finite t > \ such that supp{fi) C [—t,t], 
fi has finite moments of all orders, and dfi{x) = B{x) dx, with B a bounded non- 
negative function. Define weights w = {it'fclfe^o '"^£'^6 i'^^ on P{w). Then 
the operators F, F* , F^F , FF^ are bounded, and the operator norm for F^F and 
FF^ is bounded above by t7r|ji3||oo- 



Proof. For each k e Nq, 



Wk J-t 



x^(l){x) d^i{x) 



x''(j){x)B{x) dx 



t 

= jk J u^4>{tu)B{tu)du, 



and 



\F*^ 



^ Wk 



x''(t>{x)B{x) dx 



u''(j){tu)B(tu) du 



But by the previous lemma, this last expression is less than or equal to 

r-l 

|2 , 



So, ||-F^,||op — \\FFZ\\op 



tMBW 



= t7T\\B\ 

< tnWBW^. 



\(j)itu)\^B{tu) du 



\<f>{x)\^B{x) dx 



□ 



The following is an immediate result of Theorem 17.31 



Corollary 7.4. Suppose there exists a finite t > 1 such that supp(fi) C [— 
fj. has finite moments of all orders, and dfj.(x) = B{x) dx, with B an nonnegative 
bounded function. Define weights w = {wk} kLo where Wk '■= t^^ on ("^{w). The 
operator FF^ is an integral operator of the form 



{FF:M^)= ^_^^(y)B(y)dy. 



We can use Theorem 17.31 to study operators from the examples in Chapter [2] 
where the moments were related to the Catalan numbers. 

Corollary 7.5. Suppose /x is the Wigner semicircle measure in Examvle \2.13l 
Then under the weights w from Corollary \7.4\ FF^ is bounded, with \\FF^\\op < 47r. 
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We would also like to apply Theorem [73] to Example l2.14( but we cannot since 
the function B is not in L°°. Still, we can provide an estimate. We use a generating 
function G defined by 

(7.3) G(C) = E / 1^1' dM(x) = C'-^'^' 

where we recall the use of the notation A4k from Lemma l44l in Section |4?T1 

Lemma 7.6. Suppose supp(fi) C [— where fi has finite moments of all or- 
ders. Then the radius of convergence R of the generating function G from Equation 
(fOD satisfies R>\. 

Proof. We have 




where we use the fact that t := ess supp|a:;| on L°°(fi) and also the fact that /i is 
a probability measure. Therefore, we are assured that G converges absolutely for 
\^\t < 1. Therefore, R>j. □ 

Pick Wk := s*^ for some < s < oo. We can now establish an estimate for 
||f *(/)||^2(^) in terms of the generating function G: 



< 



|xh0(x)|d/i(a;) 



Then, using Cauchy-Schwarz, we have 



\PZ4>\\%{w) 



G\'-^Mx)\^dn{x). 



The switch of the sum and integral above is justified by a Fubini argument, since 
the power series G(x) is continuous. 

We apply this observation to Example 12.141 the secant measure. In this case, 
G{x) ~ 2(\/l — ^x"^)^^ and i = 2. Choose s > 4 and set w = {wk}^Q where 



Wk = s . Then in i (w), we have 



< 



I -if) 



M[x)\'' dpi{x) 



We next describe a result related to affine iterated function system measures. 
In particular, we show a sufficient condition under which the matrix A (see Section 
14. ip encoding an affine map r is a Hilbert-Schmidt operator. 
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Proposition 7.7. Consider the transformation t(z) = cx + h, and suppose 
|c| < 1^1 < ^- Then the operator defined by the matrix A which satisfies 

with is Hilbert- Schmidt; i.e. A is bounded, and trace{A* A) < oo. 

Proof. Let k denote the row index in A, and let j denote the column index 
in A; as before (|3.5p . we have 

^ ^({';)cHi->' o<k<j 

The («, j)*^ entry of the infinite matrix A* A is therefore 

.k) [k 



(7.4) (^*^)« = E(i) Ic'l'lb'-'l'- 



0<k<i/\j 0<k<i/\j 

and if i = j , the diagonal entries are 

k=Q ^ ' 

Set a := |cp and (i := |6p. Then, using the assumption that ^ < 1, we have 

oo j / .\ 2 



trace(AM) = J2{A*A),, (l) ^'f^'^' 



(7.5) 

oo i / .\ 2 oo ^„ • 

j=0 fe=0 ^ ^ j=Q ^ ■' 

The last series has radius of convergence j. □ 

7.2. Projection-valued measures 

This section provides an overview of the properties of projection- valued mea- 
sures, leading to the spectral decomposition for unbounded self-adjoint operators 
on a Hilbert space Ti. These ideas will be used in Section [7731 and later in Chapter [8] 
in order to discuss the Kato-Friedrichs extension of operators. For far more detailed 
treatment of this theory, we refer the reader to |Bag92| i IRS801 IRudQlj . 

Recall from Definition 11.41 that an operator H with dense domain dom(i/) in 
a Hilbert space Ti, is said to be self adjoint if H* = H and dom(iJ*) = dom(iJ). 

Definition 7.8. A Borel projection-valued measure _B on R is a map from the 
cr-algebra B of Borel subsets of M to the orthogonal projections on a Hilbert space 
H such that 

(i) Em = 0, 

(ii) if = U^iS", and S^ nS.j^9 for Si e B and i ^ j, then 



£;(5)-^i?(^,)- lini E^(^J-)' 

1=0 j=l 

where the limit is in the strong operator topology. 



70 



7. BOUNDEDNESS AND SPECTRAL PROPERTIES 



A projection-valued measure E is said to be orthogonal if 
(7.6) EiSinS2)^E{Si)EiS2) 

for all Borel sets Si, 82- We will assume all projection- valued measures are orthog- 
onal imlcss otherwise stated. 

Recall from the definition of an orthogonal projection that for every Borel set 
S, we have E{S) = E{S)* = E{S)^. 

There is a well-known theory which develops the notion of integration of mea- 
surable real- valued functions against a projection- valued measure E, resulting in an 
operator on the Hilbert space H. We will denote such an integral by f{X)E{dX). 
Quite naturally, we start by defining integrals of characteristic functions, 



E{S) = / xsWEidX), 
Jm 

then extend appropriately to measurable functions. The operator is a bounded 
operator if and only if there exists an M > such that E{\f\~^{M, 00)) = (the 
trivial projection). We say that E is a resolution of the identity if 

(7.7) In= I E{dX) = E{R). 

Example 7.9. |Bag92| Given fi a (r-finitc Borel measure on R, let H ~ Li^iy)- 
Given a Borel set S" C R, define E{S) to be a multiplication operator which multi- 
plies by the characteristic function xs'- 

E{S){J) = xsf- 

Then i? is a projection- valued measure on M, called the canonical projection-valued 
measure. Given / a measurable real- valued function on R, Jj^ /(A)£'(dA) has domain 
{.g e 7Y : fg G 7i} and on that domain, 

f{X)E{dX)(g)= fg. 



One useful property of projection- valued measures is that integration is multi- 
plicative, where the product of operators is composition. 



fiX)giX)E{dX) = j f{X)E{dX) J giX)E{dX) 

Corresponding to a projection- valued measure E, there is a parameterized family 
of real-valued measures {iyv}v<£n given by the 7i-inner product: 

(7.8) MS)^{v\E{S)v). 

Since orthogonal projections arc positive operators, the quantity {v\E(S)v) is non- 
negative for any choice of v and any set S. Given a projection-valued measure 
which is a resolution of the identity, we have 

_di.,(A) = {v\E{R)v) = \\vf. 

The following theorem gives a spectral decomposition for self-adjoint operators. 
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Theorem 7.10 ( |DS88; IRud91| ). An operator H in a Hilbert space Ti. is 
self adjoint if and only if there exists an orthogonal projection-valued measure E 
supported on the spectrum cr{H) such that 



(7.9) 



Hv 



j XE{dX) 


V = 


[ XE{dX) 






Ja(H) 



for all V G doni(iJ ) . Moreover, v G dom(iJ ) if and only if 



(7.10) 



\Hv\ 



A'^dz^t,(A) < oo, 



where the measure is defined in Equation ( |7.<^[ ). 

While we leave the details of the proof to the cited literature, it may be helpful 
here to include the computation which demonstrates the equality given in Equation 
(|7UU)) . Given v G dom(i/), 

XE{dX)v = I ( XE{dX)v [ XE{dX)v 



XE{dX)^ ^ J XEidX) 

12 \ 

A£'(dA) V V ) since the op. is self adjoint 

A^£'(dA) V v'^ multiplicative prop, of integrals 

A^di/„(A). 

□ 

The formulas from Theorem [7T0] together help us define the domains of (possi- 
bly unbounded) self-adjoint operators defined by integrals of measurable real-valued 
functions / against a projection-valued measure E. We say that is in the domain 
of an operator H = J f{X)E{dX) if 



(7.11) 



Hv- 



f{X)E{dX) 







as n — !■ oo. If H is a self- adjoint operator, we can use the equality in (|7.10p which 
gives 

r rn , 2 ^ 

X^diy^iX). 

|A|>ri 

Therefore, a given vector f G 7i is in dom(i7) if and only if 



Hv- 



XE{dX) 



(7.12) 



lim 

n — *oo 



X'^duy{X) = 0. 



7.3. Spectrum of the Kato operator 



In this section, wc show how spectral analysis of the Kato-Friedrichs operator 
translates into spectral data for the initially given moment matrix M^'') and the 
measure fi. We then illustrate the theory in several examples. Our main result 
(Theorem 17. lip is that atoms in the spectrum of the Kato-Friedrichs operator pair 
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up with atoms in the spectrum of the measure /i itself. Using this we arrive at 
our own proof of the weU- known fact (Theorem 17. 13p that the Hilbert matrix has 
continuous spectrum. 

In fact, in their spectral picture, the Kato-Friedrichs operators are directly 
related to the associated moment matrices Af*^^); this correspondence is especially 
transparent for Examples 17.191 and 17.201 In these examples, the agreement of the 
rank of the Kato-Friedrichs operators with the rank of the measure, which we define 
in Section [7. 4[ can be seen by inspection. 

Theorem 7.11. Let ji be a measure with compact support on the real line. Then 
the atoms in the spectrum of the Kato-Friedrichs operator of the moment matrix 
M = Af(^) (if any) pair up with atoms in the spectrum of the measure pL itself. 

Before we begin the proof, we will a few key points about spectral resolutions 
and operator theory. 

Let H he a. densely defined, positive self-adjoint operator on the Hilbert space 
TL. Then there exists (unique, up to unitary equivalence) a projection- valued mea- 
sure E such that 




and 

H= j \E{A\). 

J<y{H) 

In addition, for any a; G 7i, given the spectral measure Ux defined in as in Equation 
(|7.8|) . we have 

\\Hxf = j X^dvx. 

Jcr(H) 

Definition 7.12. An atom in or in £^ is a point Ai e M such that E{{\i}) ^ 
0, i.e. is not the trivial projection. 

If is a complex Hilbert space, then we can define an order on self-adjoint 
operators on 7i: 

H <K if and only if {x\Hx) < (x\Kx) for all x £ H. 

With this order, we observe that if Ei and E2 are projections on H, then 

El ^ E2 ^ El = E1E2 = E2E1. 

In this case, we say that Ei is a subprojection of E2. 

Proof. Let supp(/i) C M, where the measure fi has compact support. Denote 
the Kato-Friedrichs operator for the moment matrix Af = A'/^^^ by H. Let E be the 
projection-valued measure for H. Recall that in the unweighted case, the mapping 
F* : L^{fJ.) t"^ which takes fc{x) = X^iGNo '^^-'^^ {cijisNo is an isometry. 

Next, assume that Ai is an atom in the spectrum of H. Therefore, E[Xi) is a 
nontrivial projection, hence there exists ^1 S Range(i?(Ai)), HCilU^ = 1. We denote 
the rank-one projection onto ^1 (using Dirac notation; see Chapter [T]) by 



(7.13) 



Ex,^\^i){ii\. 
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if is a positive operator, hence Ai > 0. Since H is the Kato-Friedrichs operator, 
we have 

j %\'d^, = \\H'l^c\\l 

for aU c E v. Therefore, there exists Ai E R+ such that 

(7.14) A^Ex, < H. 

Equivalently, by the idempotent property of projections, for all c E V we have 

(7.15) A,\\ExA% < \\H'/'c\\% ^ I lApd/x. 
Recall that for all c E 

(7.16) {c\Ex,c)r-^\\Ex,c\\%^\{^,\c)e2\\ 

We now make use of the Hankcl property from Definition 11.61 in the moment 
matrix M. Let 

S : {co, ci, C2, ■■■} ^ {0, Co, ci, . . .} 
be the right shift operator in l'^. Then the Hankel condition is equivalent on V to 

(7.17) HS = S*H. 
By the spectral theorem, then 

(7.18) Ex,S^S*Ex,. 

Looking ahead to the argument in Lemma 18.81 we see that this implies that the 
vector ^1 is of the form 

ei = [l b P 63 ..f. 
Equation (|7.15p gives that for all c E T), 

2 



oo 


oo 


2 

oo 


oo 


(7.19) Ai\\Ex,c\\% ^ AiY, 








i=0 


J=0 




j=o 



Using Example 17.161 regarding the Dirac point-mass measure (5b at the point 6, 
we have 

2 

/oo 
Ifcl^dS, ^ A,\f,{b)\^ = A, Y^c.V . 

3=0 

Since the expression above is equal to the i = term from Equation (|7.19p . we also 
have 



(7.20) 



Ai J\fc\^dSh< yi/epd^ifor allcEV. 



Since the support of /i is compact, the set {/deep is dense in _L^(/.t). Therefore, 
for all Borel sets E E B{M.) there exists a sequence {cfclfegN^ C V such that 

Therefore, using the sequence {/c^} to approximate xb, we get 
(7.21) Ai I dSb< [ d^l, 
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or equivalently, 

(7.22) AiSbiE) < ^liE). 

Since Sb{{b}) = 1, it must be true that > 0. Therefore, when the Kato- 

Friedrichs operator H has an atom, the measure /i has a corresponding atom. □ 

The Kato-Friedrichs operator for the moment matrix for Lebesgue measure on 
the interval [0, 1] is the bounded operator defined relative to the standard ONB 
in £^ by the Hilbert matrix. In particular, we need not introduce weights into the 
i'^ space in order to produce a closable quadratic form. This is a case when the 
Kato-Friedrichs operator has infinite rank (see Section Fr4|) . We now sec that the 
well-known fact that the spectrum of this operator is purely continuous follows as 
an immediate corollary of Theorem lT.llI (This result was first shown in |Mag50| .) 

Theorem 7.13. The spectrum of the Hilbert matrix is continuous, i.e., there 
are no atoms in the spectral resolution of the Hilbert matrix. 

Proof. We saw that the action of the Hilbert matrix M on defines a 
bounded (see |Hal67| ) self-adjoint operator, and M is the moment matrix M*^^-* 
where fi is taken to be Lebesgue measure on the unit interval (0, 1). Moreover, the 
quadratic forms below coincide on the space D of finite vectors in i'^: 

(7.23) ll/c|U^(^) = Qm(c) = {c\Mc),2 for all c e V. 

Therefore, M is the Kato-Friedrichs operator for the quadratic form Qm- 

If M = M^f^^ contained an atom, then there would be a rank-one projection 
in the spectral resolution for M^'^\ and by Theorem 17.111 this would imply that 
Lebesgue measure restricted to (0,1) would contain an atom. That is a clear con- 
tradiction, and the therefore Af has continuous spectrum. □ 

In Section 17.51 we explicitly compute the Kato-Freidrichs operators in several 
examples and study their spectra. Recall that in general, the Kato-Friedrichs op- 
erator _ff is a self-adjoint operator on £'^{w) for some set of weights if = {wfelfegNo- 
Note that the operator H{~ H.^,) depends on the choice of w, and the choice of 
w depends on the measure ^l. The lemma below demonstrates how the weights in 
a Hilbert space £'^{w) arise in the equations governing the point spectrum of the 
Kato-Friedrichs operator. 

Recall our notation for the weighted Hilbert spaces, where lu = {w,}igNo ^re 
the weights: 

and 

{b\c)p(^) := ^ w^bid. 

In the following computations, we assume that 

(1) the measure is a positive Borel measure on M with compact support 
having moments of all orders, 

(2) the weights {wijigNo satisfy Wi > for all i e Nq, 

(3) the monomials a;* are in i^(/i) n L^{lj) for all i e Nq. 
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Lemma 7.14. Let fi be a Borel measure on M with compact support and finite 
moments of all orders. Let M ~ A/ be the moment matrix for ^, and let H be 
the Kato-Friedrichs operator for M under an appropriate choice of weights w ~ 
{wi}i^jq„. There exists X in the point spectrum of H, i.e. 

Hc^ Ac 

for some nonzero c in the domain of H , if and only if for each i CzNq, 

(7.24) ^ M,.jCj = Xw^c^. 

Proof of Lemma I7.14t The (nonnegative) real number A is in the point 
spectrum of H if and only if there is a nonzero c in the domain of H such that He ~ 
Ac. By definition of the Kato-Friedrichs operator, H is related to the quadratic form 
Q M arising from M via 

Qm{c) = ^5]c7M,-,Cfc = \\H^^^c\\l^^y Vc e V. 

J k 

When we pass to the completion, given c € dom(i/^/^) and £ L^il^), we have 
the identity 

Recall that fc G L^ifJ-) if and only if c e dom(i7^''^); this follows from the comple- 
tion in the Kato-Friedrichs construction. 

There is also a sesquilinear form Sm determined by M (connected to Qm by 
the polarization identity) which satisfies 

(7.25) 5a/(6, c) = ^ ^ 6,M,,fcCfc = (6|i/c),2(„) 

j k 

for all & e r»,c e dom(7J). 

For each i € No, let 6 = e^, the i-th standard basis vector. Equation (|7.25p 
becomes 

^MijCj = {ei\Hc)e2(.uj) = Wi{Hc)i. 
j 

We now have that He = Ac for some nonzero vector c in the domain of H if and 
only if 

MijCj = Wi{Hc)i = WiXci. 

3 

□ 

In Equations (|7.24[) . we will seek solutions for A, w, and c. We require 

c e l^{w) and ||c||^2(j„) = 1, 

i.e., 

(7.26) ^«;,|c,|2 = L 

In the examples shown in Section 17.51 we will find that it is relatively easy to 
solve (|7.26p for a single Dirac mass jj. = 5^ for any 6 G K. However, in the case of 
convex combinations of point masses, finding solutions to (|7.24p and (|7.26p is more 
tricky. 
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7.4. Rank of measures 

Consider a probability measure /i on C or on M'', and assume that fi has finite 
moments of all orders. Pick an order of the index set for the monomials, and let 
{Pk}ken^ be the associated orthogonal polynomials in L^{ijl) defined from /i; further, 
let G be the corresponding Gram matrix. Let V be the closed subspace in 
spanned by the monomials. 

Definition 7.f 5. We say that the rank ofG, and hence the rank of the measure 
fi, is the dimension of the Hilbert space V. 

We note two things in Definition l7.15l First, {PfcjfcgNg is an ONB in V. Second, 
V = L^in) if and only if the monomials are dense in L^(/x); this is known to be 
true if /X has compact support, but it need not be true in general. We will see this 
more in Chapter [S] 

Example 7.16. A Dirac measure has rank 1. 

Let 6 be a fixed complex number, and let fi :~ Sb be the corresponding Dirac 
measure. Then the moments mk of fi are 

(7.27) mfc= / z'=d/i(z) = &^ Zee No; 

Jc 

the moment matrix is 

(7.28) A^j';)=?6^ 

If V is the space of finite sequences, c = {cj}j^jqg, then 

(7.29) (c|M('')c)^2 = I J2 Cj^^ ^• 

jei% 

We showed in Example 16.61 that dim(L^(/i)) = 1, so it follows that po{z) = 1 
and pk = for k > 1. We use Equations (|4.2I1 and (|4.6p together with Lemma 14.31 
to compute the following formulas for the infinite matrices G and G^^: 



(7.30) 



G 



1 





and G"^ = 



1 
5 
62 



We see here that G has rank 1. Note that we are using the notion of inverse 
described in Section HTBl We see that GG~^ is an idcmpotent rank-one matrix. A 
direct computation further yields 



(7.31) 



G-'{G-')* ^ (b^b'") = (M^'^))*'' 



as predicted by Lemma 14.31 and Equation (|4.14p . 
Example 7.17. A measure with rank 2. 



Let n := ^{Sq + 5i). The first moment of ^ is 1, and all the other moments are 
The associated orthogonal polynomials are 

(7.32) Po{z) = l, pi(z) = 2z-l, and p2iz) ^ p^iz) = . . . = 0. 
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We then compute the Gram matrix and an inverse (in the sense of Definition 11.161 
and Example ll.lSp : 



(7.33) 



G 



1 
-12 




and G-^ 



1 
i i 

i i 



Again we can compute directly the two idempotents Ei and i?2- 

"1 



(7.34) 



and 



(7.35) 



El = GG~^ 



E2 ^ G-^G 



1 




10 

10 

10 

10 



Note that Ei is a projection; i.e., Ei = E'l = Ef, while E2 is only an idempotent; 
i.e., E2 = E2- Note that these rules are statements about infinite matrices. In fact, 
if E2 is viewed as an operator in it is unbounded (£'2(62) ^ ■) So, 

rank(G) = rank(/i) = rank(M('"') = dim(L2(^)) = 2. 



7.5. Examples 



In Example I7.18[ we use Lebesgue measure and the corresponding Hilbert ma- 
trix to illustrate the case where the moment matrix is a bounded operator. We 
also can demonstrate the IFS techniques from Chapter [3] here. Following that, we 
use the measures [i = b\ and \x = i((5o + i5i) to illustrate the spectral results from 
this chapter. In both cases, we find appropriate choices for the weights which allow 
us to find Kato-Friedrichs operators on the Hilbert space P'iui) space. Since the 
operator depends on the choice of weights, we will denote the Kato-Friedrichs op- 
erator by Hw. In Example 17. 19| for a particular choice of w, we get the associated 
Kato-Friedrichs operators to be a rank-one projection. We similarly analyze the 
spectrum of the Kato-Friedrichs operator (which has rank 2) in Example 17.201 

Example 7.18. Lebesgue measure and the Hilbert matrix. 

Let 11 be Lebesgue measure restricted to [0, 1], hence the moment matrix Af^''^ 
is the Hilbert matrix previously discussed in Example 11.71 and Section 16.11 Recall 
that it is known f |Wid66| ) that the Hilbert matrix represents a bounded operator 
on £^ with operator norm IjAff^^Hop — tt. 

It is well known that /i is exactly the equilibrium measure arising from the real 
affine IFS {Ti}i=o,i where 

■^0(2;) = ^ and Ti(a:) = 
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together with equal weights i. The IFS invariance property for moment matrices 



(7.36) Af = i (a^M^^) Ao + A* A/(^) Ai) , 

where Lemma IXTl gives the matrices and Ai (from |EST06| ): 



Ao 



1 







1 i i 

i i 

i 



Since the infinite matrices Aq and Ai arc upper triangular (and hence Aq and A^ 
arc lower triangular) . the matrix multiplication computations in (|7.36[) involve only 
finite sums. In addition, when we check the appropriate convergence issues, each 
of the matrices A^M^^f^^Ai = Af^^""^* \ for i = 0, 1, represents a bounded operator 
on e"^. We compute directly A^^Af^^Mo: 



(A*Af(^)Ao).,, = ^^(A*),.fcA4:;'(^o) 



i,3 



iAo)iA — 



W+jJ l + i 



:iAoh 



+ J 



The entries of the matrix AoA/'^'^Ao together with the IFS invariance (|7.36p prop- 
erty 



^(^A*M^f^Ao + AlM^f^Aj 



1 



1 



1 - 



1 



1 



2'+i J l + i + j 



1 



1 + 1 + j ''^ 
allow us to compute the entries of the matrix ^^Af'^^Mi: 



(AiA//(^)Ai)„- = X]^(At),,,A/(;;)(Ai) 



k=o e=o 

i j 



(7.37) 



Oi+j Z-^ 2-^ 



k=0 1=0 



1 



t \+k+l 



= 1- 



1 



1 



/ 1 + i + j 

We note that (|7.37p yields an interesting identity on binomial coefficients. Given 



i,j G No, 
(7.38) 



EE 



1 



2i+j - 1 



] 

ej i + k + £ 1 + 1 + j ■ 



fc=0 i=0 

The infinite matrices ^^^d Ai represent bounded operators on the Hilbert 
space £^ with 

\\A\\op^\\A*\Up^i i = o,i. 
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These are therefore contractive operators on 

Example 7.19. The spectrum of the Kato-Friedrichs operator for /i = Si. 

Let n = Si, the Dirac point mass at 1. We computed the Kato-Freidrichs 
operator in Example 15.151 where the weights w = {wj jjgNo C K.^ are chosen 
such that < cxD. We found that Huj is a (bounded) rank-one operator on 

£'^{w) defined by 



Recall that we also know from Example 16.61 that the dimension of L'^ifJ.) must be 
exactly 1. 

If we choose weights {wi} C such that 



E- 



= 1, 



then : P{w) £'^{w) is in fact a rank-one projection — that is, 



(7.39) 



Hi = H,, 



In this case, we know from the spectral theory of projections that the spectrum 
of is exactly {0, 1}, and the operator norm of is 1. 

If we generalize to the case where ^ = St for some real value 6, then the moment 
matrix is 

and the Kato-Friedrichs operator H^, for weights w is a rank-one operator with 
range the span of 

6 



J_ _6_ 6f_ 

Wq Wi 11)2 







Example 7.20. The Kato-Friedrichs operator associated with /i = i((5o + Si) 
When jj, = ^{Sq + Si), the moment matrix M is given by 



(7.40) 



M = - 
2 



10 





1111 
1111 
1111 



We wish to solve the eigenvalue equations (|7.24p . which in this case become 
co + ^(ci+C2H ) ^ XwqCq, fc = 



1 



(co + ci+C2H ) = XwkCk, k>l. 



(7.41) 

which can be transformed into 

1 



(7.42) 



Co ^ X{woCo - WkCk), k>l. 



by subtracting the later equations in (|7.41|) from the first. 
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To solve for c in Equation ()7.42p , we consider the cases where A = and where 
A 7^ 0. 

Case 1: If A = 0, cq must also be 0. In addition, referring back to the equations 
in the system (|7.4ip . we easily see that Equation (|7.24p is true whenever 

oo 

(7.43) ^Cfe = 0. 

k=l 

(For example, c = (0,1,-1,0,0,...).) There is no restriction on the weights 
w = {wklk&ia coming from the equations here. So for any choice of weights, 
the eigenspace for A = is infinite-dimensional. 

Case 2: If A ^ 0, the fc"^ equation for fc > in the system of equations (|7.42p can 
be solved for in terms of A, cq, and the weights w = {wklk&i'- 

(7.44) Cfe = ^U-^), fc>l. 

Wk V 2A/ 

We see that cq = implies that = for all fc, so we require cq 7^ 0. Substituting 
(|7.44p into the first equation in (|7.4ip . and cancelling cq from both sides, we obtain 
a condition on A and the weights w. 



1 °° 1 1 



Wh I \ 2A 

fc=i 



Awn. 



Define T = where this sum being finite imposes a restriction on our 

choice of weights w. Then the equation becomes 

Observe that the two conditions 

1 = Awo and 2Xwq = 1 

will lead to inconsistent systems, so wc must rule these conditions out. Otherwise, 
there are many solutions to Equations (|7.24p . 

We now demonstrate that there will be two distinct nonzero eigenvalues for the 
Kato-Friedrichs operator H^, where the weights w are chosen so that T is finite. 
Without loss of generality, take = cq = 1. Our eigenvalues A 7^ must satisfy 
Equation (j7.45p . which is now 

(7.46) '^\^i}-Y>)-^- 
This yields the quadratic equation in A: 

^"-(i + f)^ + f^»' 

This equation has distinct positive roots A+, A_ for T > 0, and they are given by 



l + ^±^/l+(? 

(7.47) A± = 

and we also see 

(7.48) A+A_ = ^ and A+ + A_ = 1 + ^. 
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Since the vector c satisfying He = Ac (where X =/= 0) is uniquely determined 
from A, Co, and w, given a choice of weights and taking cq = 1, each eigenspace 
for nonzero A is one-dimcnsional. Denote by C-|- and c_ these eigenvectors in the 
eigenspaccs for A-|- and A_ respectively. 

Going back to (|7.44p we have the simplification for the eigenvectors: 

(7.49) cfc = — =— (^), k>l. 

Wk V 2A/ Wfc V A / 

With respect to our eigenvalues X±, this gives 

/ A+-i 
c± = [I, 

then 

The reader can verify the orthogonality of the eigenvectors. 

In order to do more explicit computations, let us choose the weights w 
{uifcjfcgNo such that T = 2. This yields 

1 

which are both positive. 

Substituting back into (|7.49p we get 

(7.50) c± = (l,±^^,± 




A± =1±— , 



with 

|2 



Normalize to obtain unit vectors 



l|c±||,'2(„) =2. 



(7.51) ^± := i=c± 

which yield the rank-one projections 

E± ■■= \^±){^±\ 

onto the eigenspaces. In other words, the spectral resolution of the self-adjoint 
operator H can be written 

(7.52) H^{i + 1=)e, + {i-1=)e.. 

We next make two observations regarding the connections between H^^^, the 
square root of the Kato-Friedrichs opcrtator, and the Hilbert space i^(/i). First, 
we recall from Equation (|5.13p and Lemma 15.81 the isometry 



(7.53) y'|/,pdM=||i?^/' 



that is, the identification fc <-> H^^^c is isometric. By choosing specific weights, we 
can verify this isometry using our Kato-Friedrichs operator for fi = -^(^o + ^i). 
Let w — {wk} be given by the sequence with wq — 1 and 



Wk^2''-\ fc>l. 
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which gives T = 2. Denote F{^±) by the shorthand f±. Using ((730l) and ((73T|) . 
we get 

oo 

(7.54) /±(x)=5^(e±)fca;^ = — ± 

and, because ^ = ^(<5o + (^i), we have 



^/2 2-:e' 



1 |/±(z)|2dMx) = l±-^ 



\/2± 1 

Using Equation (|7.52p . wc have 

1/9 / 1 \i/2 / 1 n1/2 

Moreover, 

(7.55) i/i/2^±=(l±i=)'^'c± 



and hence 



V2 

which verifies the isometry in (|7.53p . 

Define the map W : L^{lj) P{w) given by 

(7.56) W{fc) ^ H^''^ciov c(^V d f{w). 

Recall from Example 16.61 that the dimension of L'^{ij) is exactly 2. Our second 
observation is that the map W is an isometry into a two-dimensional subspace of 

e{w). 

Wc will check this fact directly. Recall E± = |^±)(^±|, with ^± as in (|7.5ip . 
We only need to check that W takes an ONB in L'^in) to an orthonormal family 
in £'^{w). First, observe that the orthogonal polynomials in L'^in) arc po{x) = 1, 
Pi{x) = 2a; — 1, and Pk{x) = for k > 2. Moreover, 

Wpo = PU/(i,o,o,...) 

Ijj^ (l + -L) " 0. 0. ^ ,) + (l - ^) ' V (1. 0, 0. „) 

/ 1 \ 1/2 1 / 1x1/2 1 



and 



Wpi = W^/(-l,2,0,0,...) 



1 \ 1/2 1 / 1x1/2 1 



Since Wpo and Wpi are orthonormal in £'^{w), we have that W is an isometry. 



Example 7.21. The Kato-Friedrichs operator associated with ji = ^{Sq + 6^), 
< 6 < 1 
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This example is a more general case where ^ is a finite convex combination 
of Dirac masses, hence the associated infinite Hankel matrix will be of finite rank. 
To simplify matters we pick the atoms in /i in such a way that no weights will 
be needed. When we work with an unweighted P space, Wfc = 1 for all fc G Nq, 
and = (."^{w). The moment matrix M ~ A/'^^ for the convex combination ^ is 
equal to the Kato-Fricdrichs operator, as in the case of the Hilbcrt matrix. This 
last example also has the advantage of illustrating an ONB in L^{fi) consisting of 
rational functions; not polynomials. The choice of the rational functions is dictated 
by our Kato-Friedrichs operator H = M. 

As in the previous Example 17.201 H has a spectrum consisting of and two 
points on the positive real line, and L'^iiJ.) is two-dimensional. The point has infi- 
nite multiplicity, and the two positive eigenvalues are simple, i.e., have multiplicity 
one. 

For ji = ^{So + Sb), < b < 1, the moment matrix M is 



(7.57) 



M = H = 



1 








\b- 




¥' 


\b^ 


2" 



The orthogonal polynomials in L^{^jl) are 



(7.58) 



2x 

Po{x) = l, pi(a;) 1 -, pk{x) = 0,k>2. 

b 



However, we also have an orthogonal basis in i^(/i) consisting of the two ra- 
tional functions f±{x), where 



(7.59) 



/±(x) := a± 



bx 



1 - bx 

These correspond to the eigenvectors for H. We see that if 

^± = [a± b 62 63 
then H^± = X^± , where the parameters a± are given by 



and the two positive eigenvalues of H are 



where p is the constant 



1 




P 



1-62' 

The functions above are the images of ^± under the operator F: 

bx 



f±ix) ^ iFC±){x) = a± + 



l-bx' 
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The eigenvectors ^± are orthogonal, and it foUows that f± are orthogonal in (/x) : 

Direct computation of {f+\f-) l'^(^) also verifies that these functions are orthogonal. 

This example also demonstrates the correspondence from Theorem 17.111 be- 
tween atoms of the measure /i and atoms in the projection- valued measure E asso- 
ciated to H . The measure /i has two atoms at and 6, while the projection valued 
measure E for H has atoms at the two nonzero eigenvalues \± . 



CHAPTER 8 



The moment problem revisited 

Given a positive definite Hankel matrix M with real entries, we showed in 
Section 12.11 that there exists a measure /i (not unique in general) such that M 
is the moment matrix of ^, i.e., M = M^^^\ In this chapter, we will describe a 
setting in which the solution measure is not unique. Differing measures solving 
the same moment problem will arise from nontrivial self-adjoint extensions of a 
symmetric shift operator S. In fact, we will show that if /.j is not unique, the 
self-adjoint extensions of S yield a one-parameter family of measures satisfying the 
moment problem for M . This will yield a necessary condition for non-uniqueness of 
measure for a given Hankel matrix AI. For more details on the theory of self-adjoint 
extensions of unbounded symmetric operators and moments, see jCon90|, |RS75l 
IRud91| . 

Given a Hankel matrix M with real entries, let Qm be the quadratic form from 
Equation (|5.ip and let Hq be the Hilbert space completion of Qm- Then the inner 
product on Tig is defined on the finite sequences V by 

(cM)wq = X!5Z^»^^*J'^J■ 
^ j 

We define the shift operator S by 

(8.1) Sc = S{co, ci, ...,) = (0, Co, ci, . . .). 

The domain of S contains V which is dense in Hq, so S is densely defined. 

In Section [5TT] we observe some of the interplay between the matrix M, the shift 
operator S, and the isometry F. Then, in Sections 18.21 and 18.31 we describe how 
the self-adjoint extensions of S yield solutions to the moment problem M = M^^^. 
We also use the shift operator in Section 18.41 in order to find a Jacobi matrix 
corresponding to a given Hankel matrix M. 

8.1. The shift operator and three incarnations of symmetry 

Suppose M is a positive definite Hankel matrix and fi is a measure such that 
M = M(^). Recall that we have defined an isometry F : Hq L^{l^) by 

Fc{x) = fc{x) ^ c„x" for all ceV. 

neNo 

The shift operator S plays a major role in this chapter, and in this section we 
study how F and S behave with respect to each other. The operators F and S 
reveal three different incarnations of symmetry in the Hankel matrix M. In turn, 
these incarnations will be used in subsequent sections to prove results about non- 
uniqueness of measures which solve the moment problem, particularly in Theorem 

m 
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Incarnation 1: A symmetric operator. We see here the connection between 
the Hankel property of M and the shift operator in the Hilbert space Hq. 

Lemma 8.1. The shift operator S is symmetric in V d Tig if and only if the 
matrix M is a Hankel matrix. 



Proof. (=J>): Setting 6 = e.^, c = Cj, and {Sb\c)-HQ = {b\Sc)-n^^, we see that 



{^): Let b,ceV. Then 

{Sb\c)HQ =^^b,^iMi+jCj = ^^6,A/i+i+jCj 



(8.2) 



□ 

Incarnation 2: Multiplication by x. We notice that if p and q arc polynomials 
in L^{ii). then 



xp{x)q{x) dfi{x) 



p{x)xq{x) d/i(x). 



We can therefore state the interaction between the isometry F and the shift S. 

Lemma 8.2. Define S :~ FSF*. Then S is a symmetric operator on V d 
LP'{^), and in addition, S is a restriction (to its domain) of the multiplication 
operator : 

[M,f]{x)^xf{x). 

Proof. Let be the operator which takes f{x) to xf{x). Using the defini- 
tions of F and S, we see that 

(8.3) M^{Fc) = FSc for all ceV, 

or stated equivalently, 

xfc{x) = fsc{x) for all ceV. 

Since F is an isometry, we know that F*F is the identity on Hq. As a result, 
applying F* on the left to both sides of Equation (|8.3p yields 

S = F*M^F on V. 

Also, P = FF* is the projection of L^{n) onto the closure of the polynomials 
V C (/i) . Applying F* on the right to both sides of Equation (|8.3p yields 

M^P = FSF* ^SonV, 
which is the desired conclusion. □ 

Incarnation 3: Jacobi matrices and orthogonal polynomials. Looking 
ahead in Section [8^ we see that given the space i^(/i), there is a Jacobi matrix J 
which encodes the multiplication operator Mx in terms of a recursion relation for 
orthogonal polynomials {pk}keN„ in i^(/i). Specifically, 



J 



Po 




MxPo 


Pi 




MxPi 



8.2. SELF-ADJOINT EXTENSIONS OF A SHIFT OPERATOR 
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Lemma 18.201 then gives an intertwining relationship between the shift S and this 
Jacobi matrix J. 

8.2. Self-adjoint extensions of a shift operator 

Recall from Definition 1 1 .41 that a densely defined operator S on a Hilbert space 
H is symmetric on dom(S') C 7i if {Sh\k) = {h\Sk) for all h,k G dom(S'). Also 
recall that we can define the adjoint S* of S, as in Definition 11.31 and we call 5* a 
self-adjoint operator if S* = S**, in particular, if dom(5) = dom(5*). 

Definition 8.3. Suppose S* is a densely defined operator on I? in a Hilbert 
space Ti.. A self-adjoint extension T of 5 satisfies the following properties: 

(1) T is self-adjoint with dom(r) = dom(T*) 

(2) Tc = Sc for aU c e P. 

It may be the case that no self-adjoint extensions of S exist. 

If S is essentially self-adjoint, we call the closure of S* a trivial self-adjoint 
extension of 5*. Finally, we say that S is maximally symmetric if S has no proper 
self-adjoint extensions. 

Let S be the closure of the shift operator from Equation (|8.ip . By Lemma \8A[ 
S is symmetric in I? if and only if the matrix M is a Hankel matrix. The deficiency 
indices of S will allow us to describe its self-adjoint extensions, if any exist. 

Definition 8.4. Let S* be a closed symmetric operator on a Hilbert space H 
and let a € C with Im(a) 7^ 0. We define the deficiency subspace of S' at a by 

(8.4) £{a) ^ null(5* ~ a) = e dom(5*) : S*^ = a^}. 

Due to a beautiful argument of von Neumann, the dimension of the deficiency 
subspace will be the same for any a with Im(a) > 0, and the dimension will be the 
same for any a with Im(a) < 0. It is therefore sufficient to consider a = ±i. We 
will denote the respective deficiency subspaces for the shift operator C+ and £_ . 

Definition 8.5. The dimensions of C+ and £_ for a closed symmetric oper- 
ator S are called the deficiency indices of S and are denoted by the ordered pair 
(dim(£+), dim(£_)). Note that the indices can take on any value in No or 00. 

Definition 8.6. Let S and T be linear operators with dense domains in a 
Hilbert space, and let Gr(S'),Gr(T) be their corresponding graphs. We say that 
5'C T if Gr(S') C Gr(r). 

Note that an operator S is self-adjoint if and only if 5 C S** . If T is a self-adjoint 
extension of S, then 

S cTCT* C S*. 

Given a bounded operator J, we write JS C SJ if J maps the domain of S into 
itself and JSv = SJv for all v £ dom(5). 

Another theorem of von Neumann gives a condition under which the deficiency 
indices of a symmetric operator are equal. 

Theorem 8.7. (von Neumann, as stated in |Con901 Prop. 7.2, p. 343]) Given 
an operator S on a Hilbert space Ti., if there exists a function J : Ti. ^ Ti satisfying 
the following properties: 

(1) is the identity on Ti, 
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(2) J is conjugate linear — that is, for all a € C, J{ah) — aJ{h), 

(3) \\Jh\\ = \\h\\ for all heH, 

(4) Jdom(5) C dom(5) and JS C SJ, 

then S has equal deficiency indices, i.e. dim(£+) = dim(£_). 

The main idea in the proof of this theorem is that the operator J restricts as 
an isometry between the deficiency subspaces £+ and thus showing they have 
equal dimension. 

We apply this theorem to our shift operator S on the space Hq. Let J be the 
conjugation operator: 

(8.5) J -.Hq ^Hq with J(co,ci,C2, ...):= (c5',cr,ci', ...),. 

If ^ € then e It is readily verified that J satisfies the properties in 
Theorem 18.71 In particular we see that on V = dom(>S'), J commutes with S. 

We now calculate the deficiency indices for S, the closure of the shift operator. 

Lemma 8.8. The closed shift operator S l\8.1^ defined on D <Z H-m is either 
self-adjoint or has deficiency indices (1, 1) — i.e. 

dim{^ e dom(5*) : S*^ = i^} = l. 

In fact, if ^ & C+ with ^ 7^ 0, then iJ^ is a multiple of the vector . . .), 

where H is the self-adjoint Kato operator for the quadratic form Qm on Tig. 

Proof. If S is self-adjoint, its deficiency indices arc (0, 0). Suppose there exists 
^ 7^ such that ^ G £+. Wc will compute {H^)i via the inner product (•|-)-Hq to 
show that 

(8.6) H^eC{l,i,i^,i^,...). 

Note Ci belongs to P = dom(5) for every i e Nq. Wc now compare (-ffC)i = {^-iIOhq 
and {H^)i+i {c^+iIOhq- 

(8.7) (e.+i|0«Q = (5e,;|0«g = = *(e,;|OHQ, 
which implies by induction that is a multiple of the vector (|8.6p . 

Since H has trivial kernel in Tig, we can conclude that dim£_|_ = 1 and by 
Theorem 18. 7[ dim£_ = 1 as well. □ 

We next use a theorem, again quoted almost verbatim from |Con90| . which 
states that the self-adjoint extensions of S are determined by the partial isometrics 
from to £„ . 

Theorem 8.9. |Con90[ Theorem 2.17, p. 314] Let S be a closed symmetric 
operator. If W is a partial isometry with initial space in C+ and final space in 
then there is a closed symmetric extension Sw of S on the domain 

{f + g + Wg : / e dom(S'), 5 G initial(M^)} 

given by 

Swif + g + Wg) = Sf + tg - iWg. 
Conversely, if T is any closed symmetric extension of S , then there is a unique 
partial isometry W such that T ~ Sw o,s defined above. 

In the case of our shift operator S, we see that the only nontrivial partial 
isometrics from £+ to £_ are isometrics between the one-dimensional spaces. These 
isometrics are given by multiplication by z g C where \z\ — 1. 
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Theorem 8.10. Given a Hankel matrix M and the associated Hilbert space 
Ti.Q, let S be the closure of the symmetric shift operator jp . // S is not self- 
adjoint, then it has self-adjoint extensions Tz which have domain T) + C+ + and 
are exactly given by 

Tzic + i + Ji) = Sc + i^-izJi forzeC. 

Proof. By Lemma [8?8l if S is not self-adjoint, it has deficiency indices (1, 1). 
If ^ e then G Given z G C with \z\ = 1. we can define an isometry 
Tz : C- by = zJ^. In fact, every such isometry is of this form since the 

spaces are one-dimensionaL The resuh follows from Theorem 18.91 □ 

8.3. Self-adjoint extensions and the moment problem 

Next, we describe how the self-adjoint extensions to 5* described in Theorem 
18.101 yield solutions fiz to the moment problem M = M^'^\ (For background, see 
[ConQOj .) 

As we discussed in Section [7. 2| every self-adjoint operator T can be written in 
terms of a projection- valued measure E such that 

(8.8) T= j XE{dX). 

Js. 

PROPOSITION 8.11. [Con90( Prop. 7.2, p. 343] Given an infinite Hankel 
matrix M , suppose T is a self-adjoint nontrivial extension of the shift operator 
S on the Hilbert space TLq with corresponding projection-valued measure E ^8.8\ . 
Given E, we can define a real measure as we did in Equation \7.8^ , by 

At(-) := {eo\E{-)eo)uQ, 

where eo is the first standard basis vector (1, 0, 0, . . .). Then the real-valued measure 
jJL is a solution to the moment problem for M — that is, 

J dfi{x) = Ml for all i G No. 

Proof. From the multiplicative property of integrals against projection- valued 
measures, 

[ yE{d\) ^ T 
Jr 

for each z G No. This gives the corresponding result for /i: 

/ x^diJ.{x) = (eo|T'eo). 

JR 

Since eo G I? and T maps T> to T>, we have the following calculation: 

J x'dfiix) = (eo\T'eo)HQ = (eo|S"eo)-HQ = {eo\ei)HQ 
= Mo,j = Mi. 

□ 

We can now associate to each self-adjoint extension to the shift operator S a 
measure which satisfies the moment problem for the matrix M, in the case where 
the shift operator is not essentially self-adjoint. 
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Theorem 8.12. Given an infinite Hankel matrix M with Mo,o = 1, let S be the 
closure of the shift operator on TLq . Then the following statements are equivalent. 

(1) S is not self-adjoint. 

(2) The set of distinct solutions to the moment problem M = M^^' is a one- 
parameter family of probability measures ■ z G C, |z| = 1}. 

(3) There exist two nonequivalent probability measures fii and fi2 which are 
both solutions to the moment problem, i.e. M = M^''^^ = A/*-^^-*. 

(4) Given any measure jj. solving the moment problem M = M'^^\ the poly- 
nomials are not dense in the space 

Proof. (1) ^ (2): If S is not self-adjoint, then by Theorem 18.101 there is a one- 
parameter family {T^ : z ^ <C,\z\ = 1} of distinet (and not unitarily equivalent) 
self-adjoint extensions to S. For eaeh T^, we can define a measure /i^ as in Propo- 
sition [5T1] which satisfies the moment problem M ~ M^^^\ It remains to be shown 
that these measures are distinct. 

The isometry F : Hq L'^{ij.z) which maps T^eo to x'^ extends to map 
tp(Tz){eo) to the function tp. Thus, because the span of the functions {;^rr^ ■ ct £ 
C, Im(Q!) 7^ 0} is dense in each space Lp'{^z), the measure /i^ is uniquely determined 
by these functions. The extensions Tz are cyclic operators on TLq, which means 
(see }Con90j ) that the set of vectors {{Tz — aI)^^eQ : a G C,Im(a) ^ 0} is dense 
in Ti.Q. It follows that fiz is determined uniquely by the isomorphism F mapping 
the vector {Tz — a)~^eo £ Hq to the function jr^^ € L'^{^j). Since the operators T^ 
are distinct, this proves the measures {^z ■ z G B, |2;| = 1} are distinct. 

(2) ^ (3): Follows because each Hz is distinct. 

(3) =4> (4): We prove the contrapositive. Suppose the space of polynomials V is 
dense in £^(^) where M = M'^^ Recall the map F : T> ^ V given by Fc = J2i Ctx'' 
is an isometry on T> which extends to an isometry on Tig. Assume ^ G £+, so 
S*^ = i^ and let ceV. Then 

/ ^^F^^Fax)d^l{x) = {FSc\FOm^) 

= z(Fc|F0l2(^) =i / F^{x)FC{x)d^l{x). 
Jm 

Because F gives a one-to-one correspondence between V and V, we can say for any 
polynomial p £ V and ^ G , 

(8.10) I {x- i)^F£_{x)dfi{x) = 0. 

Jr 

The function ^ G L°°{fi), hence ^ G L^{fJ.). Let {p„} C V converge in L^(^) 
to G L'^{p). Substituting into Equation (|8.10p then gives 

/ \F£,{x)\^d^Ji{x) ^ 0, 

hence ^ = 0. Therefore, S has deficiency indices (0, 0). 

(4) ^ (1): Assume the polynomials are not dense in L^{fJ,). As before, M = M^'^\ 
Let tp G (/J.) be a nonzero bounded vector orthogonal to the subspace V spanned 
by the polynomials. Recall that L^(/x) G T' ^ if and only if V is not dense in 
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L'^{ij). Given a £ C with Im(a) ^ 0, define 



X ^ a 

Then, because each function is bounded, we have S L'^ilA- Note tliat the 
function [Mx£,a\{x) — xS,a{x) is also an L'^(^i) function because xS^a{x) — ip{^) + 

The isometry F between Hq and L^{n) resulting from Equation (|5.2p ensures 
that FF* is the projection onto the range of F, which is the closed space spanned 
by the polynomials in Therefore, we know F*'ip = 0. Given the function ^q, 

let (/)q, € Hq be defined by (pa = F*^a- We first must show that for at least one 
choice of a, 0a 7^ 0. 

Suppose that for all a e C \ R, is orthogonal to the polynomials. Define 
to be the linear span of the functions : ct G C,Im(Q:) 7^ 0}. It is well-known 

that 5* is dense in when /i is a finite real measure, so let {T/i^lfegNo C 5* be a 

sequence of functions which converges to i/j in L^(/i). This gives 



X 'ip{x)[ijj{x) - 'ipkix)]d^i{x) 

x^\ip{x)\'^d^{x) — / x^il){x)'il)kix)d^{x) 
Jr 

^ as /c 00 . 

The terms 'ijj{x)'ijjk{x) are linear combinations of functions, so the second integral 
in the sum above is zero for all k,£ G No- Therefore, for all £ e No, 

x^\ij{x)fdfi{x) = 0. 

In particular, the £ = case implies that ijj = 0. This contradicts our choice of 
ip, hence we know there must be some function which is not orthogonal to the 
polynomial space V. 
Because 



{x — a)p{x)^a{x)d^i{x) ^ / p{x)ip{x)d^i{x) ^ 
I Jr 

for all p <E "P, we have 

where Mx is the multiplication operator by x. Every polynomial p is the image of 
a finite sequence d £ T) under F. Recall that FS = M^F for all d G V, and let 
(pa = F*^a 7^ 0. Applying F* in the above equation gives 

{F*Mxp\F*^a)nQ = {Sd\(pa)nQ = a{d\(j>a)nQ ^ a{F*p\F*^a)nQ- 

Therefore, (pa satisfies the equation S*(p ~ a(p, hence (pa G C{a) for S which proves 
S is not self-adjoint. 

□ 

The following example shows a Hankel matrix M which does not have a unique 
moment problem solution. 



Example 8.13. A non-unique measure. 
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Set 

fix) / cos(xt)e-(*'+i/*'Mt. 

Then 



/.7(.)d.= (A)^ 



= for all i e No- 



Let /+ be the function max(/, 0) and let /_ be the function — min(/, 0). Set 
/ = /+-/_ and d^i±{x) f±{x) dx. Then 







.TM/i+(a;) = / xM/i_(a;). 



Next, we see that the converse of Lemma [8^ also holds. 



Theorem 8.14. Let M be an infinite Hankel matrix with Afo.o = 1- The 
following statements are equivalent. 

(1) The moment problem M = M*^^) does not have a unique solution. 

(2) Given any a e C with < Im(Q:) < 1, there exists a vector ^ G Tiq such 
that ~ A(l, a, , a^, • • • ) for some A G C. 

Proof. Let Hq be the Hilbert space completion of the quadratic form Qm 
as defined previously, and let S be the closure of the shift operator on Hq. The 
solution to the moment problem is unique if and only if the shift operator has no 
self-adjoint extensions, by Theorem 18. 121 Using Lemma [8?8l this is true if and only 
if S is self-adjoint, i.e. its deficiency indices are (0,0). 

(1 2): This is a restatement of Lemma [8.81 using a for the deficiency spaces 
instead of i. 

(2 ^ 1): Fix a G C such that < Im(Q!) < 1 and denote 

s ~ (1, a, a^,a^, . . .) G Hq. 

Assume there exists a nonzero ^ G Hq such that ~ As for some nonzero scalar 
A G C. Then for e^, a standard basis vector, 

{SeklOnc, = {Sekm)e^=Xa''+\ 

where we have selected a so that is in ^■^ as well as in Hq. We then compute 

(efelOwQ - {eklHOe^^Xa'^. 
By linearity, we have for every d £7), 

{Sd\C)nQ = OL{d\i)uQ, 

hence <^ G C{a) and we know the moment problem for ji does not have a unique 
solution. □ 

8.4. Jacobi representations of matrices 

The big picture in our work is an analysis of measures, passing from moments 
to spectra. It turns out that a number of our problems may be studied with the 
use of unbounded operators in Hilbert space, which fits in with our multi-faceted 
operator-theoretic approach to moment problems. 

In this section we focus on the special relationship between banded matrices 
which represent unbounded operators and their associated moment problems. Be- 
ginning with a Hankel matrix M, we find a (nonunique) banded Jacobi matrix T 
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which encodes the information in AI (Theorem 18. 23p . On the other hand, given a 

We 



banded matrix T, we can find an associated moment matrix M (Theorem l8.25p 
conclude with a discussion of higher-dimensional analogues and banded matrices 
which arise in quantum mechanics. 

Definition 8.15. We say that a matrix T indexed by Nq x Nq is handed if its 
nonzero entries are restricted to the main diagonal and some number of diagonal 
bands adjacent to the main diagonal. Specifically, T is banded if there exist 61,62 
such that for all j, A: S No, 

Tj,k^Q^]-bi <k<j + b2. 

When Ti and T2 are banded matrices, then 

nes 

for each pair (s, t) £ No x Nq. It is not hard to show that the banded matrices form 
an algebra with composition as the multiplication operation. 

Recall that a matrix A with real entries is called symmetric if ^ = A*'' and 
if A has complex entries, it is called hermitian ii A = A . Banded hermitian (or 
symmetric) matrices T : No x Nq ^ C define symmetric operators on ^^(No) with 
V as dense domain. 

Example 8.16. The matrix P which represents the momentum operator in 
quantum mechanics is a banded symmetric matrix: 



"0 


1 








1 













V2 





V3 



.11) 



p 



V^n~T 





Given v & i'^, we have 







The following definition generalizes to higher dimensions, but for clarity we 
will state everything in a one-dimensional form. We begin with a Hankel matrix 
M, and again work in the Hilbert space Hq. We denote the standard orthonormal 
basis in ^^(Nq) by {(5n}neNo- Even though {5„}„gNo is not an ONB in Hq, we also 
use Sn to denote the element of Hq with O's in every place except the {n + I)''*, 
which contains a 1. 

Definition 8.17. We say that a Hankel matrix M is of T-class if there exists 
some banded hermitian matrix T, representing a symmetric operator, such that 

(8.12) Af,-fc = {5^\T^+H^)p. 
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Recall from Proposition lS. 1 ll that if M is Hankel and positive semidefinite, then 
M satisfies the equation 

where is a self-adjoint extension of the shift operator on the Hilbert space TLq. 
We will demonstrate the correspondence between this shift operator S (when a 
moment matrix M is given) and a symmetric Jacobi matrix T. 

Lemma 8.18. Every Hankel matrix of T -class is positive semidefinite. 

Proof. Let v £ V. Then because T is banded, all summations are finite: 



.13) ={J2v,T^6o\J2v,T% 

j k 



□ 



Lemma 8.19 (One-dimensional version). Suppose M is Hankel of T-class. De- 
fineW -.HQ^e (No) by 

W{5k) := T'^So Vfc e No. 

Then W is an isometry. 

Proof. We first show that |jW(<5fe)|j^2 = ll<^fellwQ- Recall from the norm on 
Hq that 



4llwo = Qa/(4), 



and 

QM{6k)^Mk,k = {So\T''+''6o)i2. 

On the other hand, 

\\W{5k)\\% = {tHo\T^5o)p- = {5q\T^+Ho)p. 

Since every element w of 2? is a finite linear combination of 8k&, the rest of the 
proof follows from the same computation used in the proof of Lemma 18.181 □ 

As in Equation (|8.1[) . let S will refer to the closure of the shift operator, which 
is defined on V by 

S{co, ci, C2, ...):= (0, Co, ci, C2, . . .). 

Because we will consider deficiency indices of two different operators S and 
T, we let C±{S) denote the dimension of the deficiency subspace of 5", and we let 
C±{T) denote the dimension of the deficiency subspace of T. 

Lemma 8.20. Let M be of T-class and let the Hilbert space Hq be as defined 
above. On the dense domain T> C Tig, WS = TW . 
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Table 1. Relationships among S, T, and W. 
S r\ r\T 

Proof. Let c = (co, ci, . . .) e P. Then 

WS{c)^W{0,co,c,,...) 

.14) =coTSo + ciT^So + C2T^So + --- 

= T{coSo + ciTSo + C2T^So + ■■■ = TW{c). 

An immediate resuh of Lemma 18.201 is the following: 



□ 



Lemma 8.21. Suppose the set {Wc}ceT> *s dense in dom{T). Then the isometry 
W maps the S -defect subspaces C±{S) into the T -defect subspaces C±(T). 

Proof. Suppose f± denotes an element of C±{S). We know that S* f± = ±i/± 
if and only if {Sc\f±)'HQ ~ zki{c\f±)'HQ for all c E V. With the assumption of 
density, T*Wf± = ±iWf± if and only if 

{TWc\Wf±)i2 = ±i{Wc\Wf±)i2 

for all c e JF. 

Let c e P. Then by Lemmas KM and KW\ 

{TWc\Wf+)i2 = {WSc\Wf±)e2 = {Sc\f±)nQ 
^^■^^^ = {c\S*fi)nc, - ±*(c|/±>«« = ±t{Wc\Wfi),2. 

□ 

Recall that we established in Lemma 18.81 that the deficiency indices of S are 
either (0,0) or (1,1). 

Corollary 8.22. Suppose the set {Wc \ c G J-^} is dense in dom[T). Then 

C±{T)>C±{S). 

Table [1] shows the relationships among S, T, and W. 
Because 

(8.16) TW = WS, 

when S{co, ci, . . .) := (0, cq, ci, . . .) in Timin, we can take adjoints in (|8.16p to see 

that 

(8.17) W*T* = S*W*. 

Finally, W* maps C±{T) into C±{S). To see this, suppose T*/± = ±«/±- Apply 
W* and use ([STT]) : 

W*T*f±^±iW*f±= S*W*f±, 

which implies that W* f± G C±{S). 

We now come to the main result of this section: every positive definite Hankel 
matrix M is of T-class, and the matrix T such that 

M,. fc - {6o\T^+''6o)i2 
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can be chosen to be banded with respect to the canonical ONB in £^(No). 

Theorem 8.23. Let M he a positive definite Hankel matrix, where M is nor- 
malized so that Afo,o ~ 1- 

(a) Then there exists a banded symmetric matrix T operating on P(Nq) such that 
118. is satisfied. 

(b) We may choose T of the handed form 

'bo ao ■■■ ' 

oq bi ai 

oT &2 ^2 

ai &3 as 



(8.18) 



T : 



Proof. By Theorem 12.21 there exists ^ such that 

JR 

Now select the orthogonal polynomials po{x),pi{x),p2{x), ... in L'^{pi) yjiih.po{x) 
1. Then the mapping 



.19) 



fc=0 



is an isometry of a subspace in L^(^) onto £^(No). Specifically, 

2 



.20) 



k k 

It is well-known (see, for example |Akh65l IAAR99] ^ that there exist sequences 
{flnlneNo s-nd {6n}No such that the following three-term recursion formulas are 
satisfied: 

xpo = fooPo + aopi 

xpi = aopo + bipi + arp2 



xpj = Qj-ipj-i + bjPj + ajpj+i 



Because of the Hankel property assumed for M, in proving (|8.12p . it is enough to 
show that 

(8.21) Mo,fc = ((5o|r%),2. 

Using ((51^ and ^^W^, we have for all fc G Nq, 

{So\T''So)i!2 = / pox''podfi{x) 



x'' d^i{x) = Mo,fc, 



which is the desired conclusion. 



□ 
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Remark 8.24. There are many other choices of banded symmetric or hermitian 
matrices which solve (|8.12p . The candidates for T are dictated by applications. 

Consider a fixed positive definite Hankel matrix M such that the associated 
symmetric operator 5* has deficiency indices (f , f ). Then by Theorem lS.f 2[ there is a 
one-parameter family of inequivalent measures : \z\ = f } such that M^'^-' = M. 
Further, for each measure /x^, we may compute an associated symmetric Jacobi 
matrix , as in Theorem 18.231 The question we ask next is whether the Jacobi 
matrix depends on the measure /iz used to compute the orthogonal polynomials. 

Consider the Hilbert space Tig, in which the shift operator S has dense domain. 
The Jacobi matrix we just computed is also a symmetric operator, this time with 
dense domain in ^■^(Nq). Moreover, for each /i solving the Af-moment problem, we 
have an isometry F ~ which maps Tig into (/i) and which intertwines S with 
multiplication by x ('Lemma l8.2p . Recall also that encodes multiplication by x. 

By looking at the relevant formulas for the orthogonal polynomials {pn} C 
LP'{^z), we see that does not depend on z. In particular, we use the following 
well-known formulas for the orthogonal polynomials to justify our answer [Akh65j . 

Define 

mo ... mfc 
mi ... TOfc+i 



Dk = det 



mfc 



and let Sm = {/s G N : Dfe 7^ 0}. Then if A: e 5*^/ 
where 

mo 



Dk{x) = det 



mi 



rrik-i 
I 



mi 
m2 



mfc 

X 



m2k 

then pk{x) = (Dfe- 



rrik+i 



m2fc-i 



,Dk)-'/^Dk{x), 



As a result, the isometry is also independent of ^; the closed subspace 
spanned by the polynomials in LF'{^,z) is independent of z, and it is only the relative 
orthogonal complement in L^(/iz) which depends on z. The orthogonal complement 
can be described by 

It/) G L^iix) j iix)x''dix{x) = for aU fc G No|. 

Lemma 18.21 tells us that S and T are unitarily equivalent. Moreover, 

so the spectral measures derived from the self-adjoint extensions S oi S and T of 
T produce the same measures /i which solve the M-moment problem. Thus, if 
we are thus given the Jacobi matrix T, the self-adjoint extensions of T in ^^(Nq) 
correspond to spectral measures /i which will have a common moment matrix M . 
We summarize with a theorem: 

Theorem 8.25. Let M be a positive definite Hankel matrix, and let fi be a 
measure such that M^'^) = AI . The associated Jacobi matrix T is independent of 
the choice of measure n which solves the M-moment problem. Conversely, every 
symmetric Jacobi matrix T gives rise to a moment problem. 
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The following three conditions are equivalent: 

(a) The solution fi to the M -moment problem is unique. 

(b) The deficiency indices for T are (0,0). 

(c) The polynomials are dense in 

Proof. The forward direction is shown in the discussion preceding the state- 
ment of the theorem. Given a symmetric Jacobi matrix T, define a Hankel matrix 
M by taking 

Since powers of banded matrices remain banded, the entries of M are all finite. 
By Lemma 18.181 M is a positive-definite Hankel matrix, hence yielding a moment 
problem M = Af '^'^^ which has at least one solution. 

The equivalent statements follow from Theorem 18.121 Corollary [8221 and the 
discussion above showing that C±{T) < C±{S). □ 

8.5. The triple recursion relation and extensions to higher dimensions 

Let poix) = l,pi{x),p2{x), ... be the orthogonal polynomials with respect to 
p., where deg{pk) = k. (Here we assume that /i corresponds to a positive definite, 
not positive semidefinite, linear functional on L'^{p).) We can use Gram-Schmidt 
on {1, x,x^ , . . .}, so that 

span{l,a;,a;^, . . .} = span{po(a;),Pi(a;),P2(a;), ■ • ■} 

and 

j pj{x)pk{x) dp{x) ~ Sj^k (the Kronccker delta). 
Then for fc G No we have 

(8.22) xpn~i[x) = &,i-lP,i-2(a;) + an~iPn-i{x) + bnPn{x) 
and 

(8.23) Xpn{x) = bnPn~l{x) + OnPnix) + 6„+ip„+i (x). 

Note that the 6„ appearing in (|8.22p is the conjugate of the 6„ appearing in (|8.23p . 
To see this that this is true, suppose 

xpnix) = Bpn^i{x) + terms in p„(x) and p„+i(x) 

xpn-i{x) = Cpn{x) + terms in p„-2{x) and Pn-iix). 

We now take advantage of the orthogonality: 

B = {pn~l{x)\xpn{x)) L2^^) = {xpn-l{x)\pn{x)) L^^,) 

= (Cp„(a;)|p„(a;))i2(^) = C(p„(a;)|p„(x-))i2(^,) = C. 

The triple recursion relation can be considered in terms of projections onto finite- 
dimensional subspaces. Set 

Hn := span{l,x, . . .^x"} = spa.n{po{x),pi{x), . . . ,Pn{x)} 

and let Qn be the projection onto 7i„ with 

Qn Qn Qn 

and 

Qn{L\fl))=nn. 
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(Note: Qn is a projection, and is not related to our earlier quadratic form.) We 
are interested in the shift operator, which is the same as multiplication by x on 
V C L^{^J.)■ We have 

xHn C Hn+l, 

so 

Set Q;^ := I — Qn — that is, is the projection onto 7i^, where Ji^ — (fi) QTLn- 
If A: < n — 1, then (p^-, xpn)L^{fj,} — 0. To see this, write 

and pn is orthogonal to any polynomial with degree less than n. 
Now, 



xp„ 



i{xpn) + Qni^Pn) 



hnPn~l + CLnPn + bn+lPn+1 



If we restrict J* to V, the subspace of all polynomials, then 
since J C J* . Moreover, since the projections Q„ are self-adjoint, 

Qn + lJQn — JQn QnJQn-\-l — QnJ • 



Finally, 



since 



Qn~2{JQn) — Qn-2JQn-lQn — Qn-2JQn~l- 



The projection approach to the recursion relation can be extended to M'' for 
d > 1. For each of the d coordinate directions, there is a Jacobi matrix 



Jk = 



"o 

















bi'^ 











"2 















Ak) 

"3 


[k) 


Ak) 

04 













(fc) 



and a shift operator 5^., which is realized by multiplication in the fc*'' coordinate 

Skp{xi, ...,Xd)= Xkp{xi, ...,Xd). 

Recall that the degree of x" = x^^ ■ ■ ■ x'^'' is ai + . . . + ad- With this notation, the 
finite-dimensional subspaces are 

Hn = {pe V\degip) < n}. 
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Table 2. Two approaches to moments and banded matrices. 



Banded T 


— > 


Hankel Mr 


— > 


measures /i 



Hankel M 


— > 




— > 


Banded Tm 



8.6. Concrete Jacobi matrices 

We examine the momentum and position operators from quantum mechanics 
for one degree of freedom and then study some Hamihonian operators (for exam- 
ple, the polynomials in the momentum and position operators in Table [3]). The 
operators we consider all have a common property: in a natural orthonormal basis 
their representations take the form of infinite banded matrices; that is, the matrices 
have zeros outside a band around the diagonal of finite width. 

The banded property of the matrices makes matrix multiplication easy; under 
multiplication, the banded matrices form an algebra of unbounded operators. While 
such banded matrices follow simple algebraic rules, their spectral theory can be 
subtle. For example, we show that these operators may not have a well-defined 
spectral resolution. Using von Neumann's deficiency indices, we showed above the 
connection of Jacobi matrices to the theory of extensions of symmetric operators 
with dense domain, and thereby to moment problems. 

One often encounters problems in physics, such as the Heisenberg banded ma- 
trices T, where the nature of the bands is dictated by applications. A particular 
infinite matrix T represents an operator in an P sequence space. In fact, in a par- 
ticular application, T may be realized in a different Hilbert space, for example an 
function space, but the function version will be unitarily equivalent to the ma- 
trix model. In particular, we refer to the fact that Heisenberg's matrix formulation 
of quantum mechanics is unitarily equivalent to Schrodinger's wave formulation 
in function space. For example, the momentum operator P is represented by the 
matrix (|8.1ip in P and by the operator i-^ on a dense subspace of L^(IR). See 
Examples 18. 161 and 18.261 and the remark following the two examples. 

From our banded matrix T we then get a Hankel matrix A/, and we apply 
our theory to M. In particular, we find the measures ^ which solve the moment 
problem for M . We use operator theory in constructing the family of measures /i 
which solve the moment problem at hand. 

We emphasize further that in applications, one typically encounters a much 
richer family of banded symmetric or hermitian infinite matrices T; for example 
those from Heisenberg's quantum mechanics. In these matrices, the band-size will 
typically be more than three. In fact the band can be any size, and the deficiency 
indices can be anything. However this wider class of banded matrices, including 
for example anharmonic oscillators, may be studied with the aid of the associated 
Jacobi matrices. 

We begin with two Jacobi matrices, the matrix P given in Example 18.161 and 
Q here. 

Example 8.26. The operator Q is represented by multiplication by x on L^(]R). 
The operator Q can be represented on £^ by a matrix defined by 

{Qv)n = - lVn-1 - \/nVn+l)- 
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Remark 8.27. The two operators P and Q in Examples 18.161 and 18.261 have 
dense domains in L^(]R): 

dom(P) = {/ e L'-{R) : /' G L^R)} 

and 

dom(Q) = {/ e L^iR) : xf{x) e L^{R)}. 
The operators P and Q are both self-adjoint, and they are unitarily equivalent via 
the Fourier transform in L^(R). 
Setting 

A± := -^{P±iQ), 
we get A*^ = , and the commutator 
(8.24) [A+,A_] = -/. 

The Hermite function /io(a;) = cqb^^ satisfies 

A^ha = 0. 
An application of Equation (|8.24p yields 

/lo = nA"/io for each n e N. 
The functions /i„ ;= c„A+/io diagonalize the harmonic oscillator Hamiltonian 

H ■.= A+A^ = ^{P'+Q'' -I), 

and the constants c„ can be chosen so that {hn\n e Nq} is an orthonormal basis in 
L^(R) consisting of Hermite functions. 

Using this ONB, we arrive at the two matrix representations for P and Q in 
Examples 18.161 and 18.261 Specifically. 

{hn-l\Phn) = ^\/n, 
(h,,\Phn) =0, 

and ^ 

{hn+l\Phn) = -\/n+ 1, 

which is the Jacobi matrix in (|8.1ip . 

Example 8.28. If T = QPQ in ^^(No), then T has deficiency indices both 
equal to 1. 

Proof. A differential equations problem. □ 
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Table 3. Some polynomials in the position and momentum operators. 



T 


index 


spectrum 


QPQ 


(1,1) 


depends on the choice 
of selfadjoint extension 




(0.0) 


discrete, anharmonic oscillator 


P'-Q' 


(2,2) 


repulsive potential 
quantum particle shoots to infinity in finite time 


P^ + Q^ 


(0,0) 


{2n+l 1 71 e No} 


P^~Q^ 


(0,0) 


continuous M (easier to see in _L^(M) 
than by using matrix calculations) 
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